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Introduction 


The  constant  development  in  wireless  system  technologies  enables  inovative  new  applications 
that  introduce  new  challenges  and  demands  on  antenna  engineers.  These  trendlines  show  need 
for  advanced  smart  antenna  systems  in  communications  and  surveillance  with  a  variety  of 
beamforming  and  beamsteering  capabilities.  Additionally,  electromagnetic  performances  of 
antenna  have  often  to  comply  with  physical  system  requirements  such  as  customized  shape 
and/or  lightweight.  Curved  conformal  microstrip  arrays  have  a  potential  to  meet  all  these 
requirements. 

The  antenna  arrays  that  conform  to  cylindrical  structure  can  enable  radiation  pattern  with 
directed  beams  in  arbitrary  directions,  or  omnidirectional  pattern  in  plane  perpendicular  to  the 
cylinder.  With  their  cylindrically  shaped  structure  they  can  be  easily  mounted  on  various 
airborne  objects  (fighter  jets,  smart  missiles,  rockets)  without  disturbing  their  aerodynamical 
properties.  On  the  other  hand,  spherical  arrays  can  steer  single  or  multiple  beams  throughout  the 
complete  hemisphere  without  loss  of  gain  in  any  direction.  Therefore,  spherical  array  is  good 
candidate  for  satellite  terminals,  telemetry  and  command  applications. 

As  long  as  the  radius  of  the  curved  structure  is  large  compared  to  the  wavelength,  the  antenna 
can  be  approximated  with  the  planar  structure  of  infinite  extent.  However,  this  is  not  the  case  for 
the  most  of  the  practical  conformal  antennas.  Therefore,  computational  tools  are  required  which 
enable  fast  and  accurate  analysis  of  modern  conformal  air  interfaces.  Additionally,  in  case  of 
array  applications,  it  is  critical  to  be  able  to  take  into  account  the  mutual  coupling  between  the 
array  elements.  In  more  detail,  the  mutual  coupling  can  cause  significant  change  in  element 
pattern  of  array  elements,  and  therefore  must  be  included  in  the  design. 

In  course  of  this  project  a  computer  program  is  developed  that  performs  a  rigorous  analysis  of 
spherical  and  cylindrical  arrays  of  microstrip  patch  elements.  The  dielectric  structure  of  the 
antenna  or  the  antenna  array  is  rigorously  taken  into  account  by  using  proper  Green's  functions 
in  spectral  domain.  The  Method  of  moments  (MoM)  is  applied  to  analyze  array  by  means  of 
element-by-element  approach,  thus  taking  the  mutual  coupling  into  account  as  well.  All 
important  information  about  the  antenna  array  can  be  calculated  by  the  software:  the  input  port 
impedance,  radiation  pattern  as  well  as  the  mutual  coupling  between  the  elements  can  be 
obtained. 
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The  capabilities  and  accuracy  of  the  computer  program  are  tested  on  different  practical  designs. 
One  of  them,  the  spherical  test  model,  is  developed  for  this  purpose.  The  comparison  of  some 
spherical  and  cylindrical  designs  with  their  planar  counterparts  is  given  in  order  to  show  impact 
of  curvature  on  input  impedance  and  radiation  properties.  All  of  these  results  will  be  presented  in 
this  report. 


2  PROBLEM  IDENTIFICATION 
AND  ITS  SIGNIFICANCE 
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Problem  Identification  and  its  Significance 


Conformal  antennas  are  defined  as  antennas  that  conform  to  the  structure  for  other  reason  than 
electromagnetic.  This  reason  can  be  aerodynamic,  hydrodynamic,  aesthetic  or  coverage  reason. 
Because  of  these  potential  advantages  over  planar  antennas  conformal  microstrip  patch  antennas 
have  been  proposed  for  a  wide  range  of  applications.  For  example,  future  military  airborne  and 
space-based  sensors  will  require  active  electronically  scanned  array  antennas.  Planar  arrays  are 
possible  solution.  However,  planar  arrays  have  inherent  disadvantage  because  they  can  be 
electronically  scanned  to  only  about  60  degrees  fi^om  boresight,  and  the  gain  falls  with  scanning 
angle.  The  conformal  arrays  can  scan  the  beam  without  drop  of  the  gain,  and  the  size  of  the 
conformal  antenna  needed  for  obtaining  the  same  gain  as  its  planar  counterpart  is  approximately 
the  same.  Spherical  arrays  have  possibility  of  directing  single  or  multiple  beams  through 
complete  hemisphere  and,  therefore,  spherical  arrays  are  an  attractive  solution  for  satellite 
tracking,  telemetry  and  command  applications.  Other  example  is  cylindrical  microstrip  arrays 
mounted  on  smart  missiles  and  rockets  without  disturbing  aerodynamic  properties. 

Conventional  phased  array  antennas  are  thick,  heavy  and  expensive,  and  require  large  amounts 
of  power.  Future  military  airbone  and  space-based  phased  arrays  will  use  conformal,  low-weight 
and  low-power  elements  in  order  to  meet  extremely  stringent  system  requirements.  Therefore, 
new  solutions  are  proposed  for  such  applications.  One  of  the  most  promising  is  to  use  Micro 
Electro  Mechanical  Switches  (MEMS)  as  an  essential  part  of  the  low-loss  phase  shifters.  The 
proposed  architecture  for  the  building  block  of  phased  arrays  is  a  sandwich  structure  where 
microstrip  antenna,  RF-distribution  layer,  T/R-switch,  low  noise  amplifiers,  power  amplifiers 
and  phase  shifters  are  built  in  a  microstrip  multilayer  structure.  The  radiating  element  in  this 
geometry  is  an  aperture-coupled  microstrip  patch  antenna. 

The  traditional  approach  of  analyzing  conformal  microstrip  antennas  is  to  approximate  the 
conformal  structure  with  locally  planar  one,  and  then  to  use  some  method  for  planar  antennas. 
This  is  a  reasonable  approximation  if  the  radius  of  the  structure  curvature  is  large.  However,  for 
small  radiuses  the  properties  of  the  antenna  begin  to  differ  significantly  from  their  planar 
counterpart. 

Some  differences  between  planar  and  cylindrical  patch  antennas  are  as  follows.  For  axially 
polarized  patch  antennas  the  input  impedance  will  significantly  drop  with  decrease  of  the 
structure  radius.  For  patches  polarized  in  circumferential  direction  the  change  of  input 
impedance  is  opposite,  i.e.  input  impedance  will  grow  with  decrease  of  the  structure  radius.  The 
amplitude  of  surface  waves,  and  therefore  mutual  coupling,  is  enlarged  in  axial  direction  and 
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reduced  in  circumferential  direction  in  comparison  to  the  planar  case.  The  radiation  pattern  can 
be  significantly  changed  in  circumferential  direction.  For  example,  we  can  get  omnidirectional 
radiation  pattern  if  the  patch  is  almost  wrapped-around,  where  in  the  planar  case  the  backward 
radiation  is  always  weak.  The  resonant  frequency  is  almost  the  same  for  axially  polarized 
patches,  and  it  is  significantly  changed  for  patches  polarized  in  circumferential  direction. 

Spherical  antennas  have  different  properties.  The  resonant  frequency  is  changed  a  lot  with 
decrease  of  the  structure  radius,  while  the  resonant  resistance  stays  almost  constant.  The  mutual 
coupling  effects  are  reduced  in  comparison  to  the  planar  case.  Like  in  the  cylindrical  case,  the 
radiation  pattern  can  be  significantly  modified.  As  expected,  the  back-radiation  is  smaller  for 
spheres  with  larger  radius. 

The  significant  differences  in  properties  of  conformal  antennas  show  that  programs  for  analyzing 
planar  microstrip  antennas  are  not  sufficiently  accurate  for  design  of  conformal  antennas,  and 
thus  there  is  a  need  for  programs  for  analyzing  conformal  microstrip  antennas.  Such  programs 
will  rigorously  take  into  account  the  effects  of  the  curvature,  and  therefore  will  help  the 
engineers  in  designing  the  conformal  antennas. 

Uniqueness:  It  was  planned  to  develop  two  programs  as  outcomes  of  the  project:  (a)  program 
for  analyzing  patch  arrays  printed  on  spherical  structures,  and  (b)  program  for  analyzing 
aperture-coupled  patch  arrays  on  cylindrical  structures  (the  upgrade  of  the  existing  program 
“CyMPA”).  Spherical  patch  antennas  are  rarely  considered  in  scientific  literature.  Till  now  only 
a  single  patch  of  circular  or  ring  shape,  fed  by  a  coax-line,  on  single-layer  spherical  structure  is 
analyzed  (state  of  art  about  this  subject  can  be  found  e.g.  in  recently  published  book:  K.L.Wong, 
"Design  of  Nonplanar  Microstrip  Antennas  and  Transmission  Lines”,  Wiley  1999.).  Therefore, 
the  development  of  the  program  for  analyzing  arrays  of  rectangular  or  circular  patches,  where 
mutual  coupling  between  elements  are  rigorously  taken  into  account,  is  scientifically  new.  The 
chosen  analysis  method  for  aperture-coupled  conformal  patches  is  extension  of  the  method 
which  is  used  for  planar  geometries,  and  which  in  planar  case  shows  good  agreement  between 
calculated  and  experimental  results. 
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3  PROJECT  OBJECTIVE 
AND  REALIZED  OUTCOMES 
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Project  objective  and  realized  outcomes 


We  proposed  a  two-phase,  24-month  effort  to  develop  a  model  for  analyzing  single-curved  and 
double-curved  conformal  microstrip  antennas.  It  was  planned  to  develop  two  computer  programs 
as  outcomes.  First  proposed  program  is  for  analyzing  patch  arrays  printed  on  spherical 
structures,  where  the  patches  are  of  rectangular  or  circular  shape  and  the  feeding  structure  is 
microstrip  transmission  line,  coaxial  transmission  line,  or  aperture-coupled  feeding.  The  interest 
for  such  a  program  is  due  to  the  fact  that  conformal  arrays  in  some  applications  cannot  be 
approximated  as  cylindrical  ones,  i.e.  as  structures  with  one  curvature  direction,  and  thus  another 
curvature  direction  should  be  taken  into  account.  Second  proposed  program  is  an  upgrade  of  the 
program  "CyMPA"  that  analyzes  patch  arrays  on  cylindrical  structures  (i.e.  on  single-curved 
structures).  The  upgrade  will  include  the  analysis  of  aperture-coupled  patch  arrays  on  cylindrical 
structures. 

It  was  planned  that  both  programs  calculate  the  following  antenna  characteristics: 

•  current  distribution  at  each  patch  in  the  array 

•  input  impedance  at  each  input  port  in  the  array 

•  mutual  coupling  between  each  two  patches  in  the  array 

•  radiation  pattern  of  the  array  when  all  mutual  couplings  are  taken  into  account 

•  radiation  pattern  of  the  array  without  taking  mutual  coupling  into  account  (fast 
calculations  of  the  radiation  pattern  needed  for  making  first  design  of  the  array). 


During  first  12  months  phase  we  have  developed: 

(a)  Program  “SMiPA”  that  analyzes  patch  arrays  printed  on  spherical  structures.  The  patches  are 
of  rectangular  or  circular  shape,  and  they  can  be  fed  by  microstrip  or  coaxial  transmission 
line.  The  program  calculates: 

•  current  distribution  at  each  patch  in  the  array 

•  input  impedance  at  each  input  port  in  the  array 

•  radiation  pattern  of  the  array  without  taking  mutual  coupling  into  account. 


Analysis  of  Curved  Microstrip  Antennas 


11 


(b)  The  upgrade  of  the  program  “CyMPA”  that  calculates  aperture-coupled  patch  antennas  on 
cylindrical  structures.  Therefore,  with  the  upgraded  program  it  is  possible  to  fully  analyze 
aperture-coupled  cylindrical  patch  arrays  (input  impedance,  mutual  coupling,  radiation 
pattern  with/without  considering  mutual  coupling). 


The  results  of  the  second  12  months  phase  are: 

(a)  We  have  finished  the  program  for  analyzing  spherical  patch  arrays.  The  newly  developed 
part  of  the  program  calculates; 

•  mutual  coupling  between  each  two  patches  in  the  array 

•  radiation  pattern  of  the  array  when  all  mutual  couplings  are  taken  into  account 

•  aperture  coupling  as  a  feeding  possibility. 

(b)  We  have  developed  an  experimental  model  for  testing  the  program  for  analyzing  spherical 
arrays  since  there  are  only  a  few  experimental  results  in  literature  in  connection  with 
spherical  patch  antennas. 


The  developed  programs  are  written  in  FORTRAN  program  language  since  there  is  no  faster 
program  language  for  computational  physics  (computational  electromagnetics).  Program  is 
independent  on  the  machine,  and  it  can  be  compiled  and  run  on  every  machine  that  has 
FORTRAN  90  compiler.  Furthermore,  for  MS  Windows  PC  environment,  a  graphical  user 
interface  is  developed  to  make  setting  of  input  parameters  easier,  and  to  obtain  graphical 
presentation  of  the  results.  On  other  machines  or  operating  systems  communication  with  the 
program  is  made  via  input/output  ASCII  files.  In  more  details,  the  input  file  should  be  filled 
before  running  the  program,  and  the  results  are  written  into  the  output  file  that  can  be  graphically 
presented  by  any  graphical  program. 


4  PROJECT  OUTCOMES 
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4.1  Analysis  of  microstrip  patch  antennas 
on  spherical  structures 


4.1.1  Introduction 

The  purpose  of  this  chapter  is  to  describe  a  program  that  analyzes  microstrip  patch  arrays  on 
spherical  structures.  The  patches  can  be  of  rectangular  or  circular  shape,  and  the  spherical 
dielectric  structure  can  be  multilayered.  The  program  calculates  radiation  pattern  of  the  array, 
input  impedance  of  each  patch  element  and  mutual  coupling  coefficients  between  each  two  patch 
elements.  The  patches  can  be  fed  by  coaxial  or  microstrip  transmission  line,  or  they  can  be 
aperture-coupled  (the  model  for  analyzing  aperture-coupled  patches  is  described  in  next 
chapter).  The  structure  is  rigorously  taken  into  account  by  using  proper  Green’s  functions,  and 
the  electric  field  integral  equation  is  numerically  solved  by  applying  moment  method. 


4.1.2  Method  of  Analysis 

The  geometry  of  the  problem  is  shown  in  Figure  4. 1 .  The  rectangular  or  circular  patches  are 
placed  on  or  embedded  in  a  multilayer  spherical  structure.  The  radius  of  the  grounded  sphere  and 
of  the  patches  are  Vond  and  rpatch,  respectively.  The  dimensions  of  each  quasi-rectangular  patch 
are  fVe  and  (fVe  =  2  6^,  r patch ,  =  2^pr patch)  and  the  dimension  of  each  circular  patch  is  We 

(We  =  2  Qprpatch)-  The  distances  of  the  feeding  point  from  the  center  of  the  patch  in  the  0  and  (|) 
directions  are  d/eed  and  ^feed ,  respectively.  Two  types  of  feeding  structures  are  analyzed  in  this 
chapter:  coax  line  feeding  and  microstrip  line  feeding. 

To  determine  the  current  on  the  patches,  we  consider  the  integral  equation  for  electric  field 
components  tangential  to  the  patches  [1]: 


(  E)nc(«I/fecrf)  ^scal(^patch))  taii~0  (1) 

Here  J/eed  and  Jpaich  denote  the  current  on  the  patch  and  the  feeding  structure,  respectively.  The 
unknown  patch  current  is  expanded  into  sum  of  the  basis  functions: 

N 

^ patch  ~  (2) 

/=! 

The  unknown  coefficients  a,  are  determined  by  applying  the  moment  method  (MoM).  We  have 
used  the  same  test  functions  as  basis  functions  (Galerkin’s  method). 


(3) 


N 

~  (*J  J  > ^scat  (.^i ))  ~  ’ ^inc  feed ))  • 

/=1 

The  elements  of  the  impedance  matrix  [Zji\  and  voltage  vector  [Vi]  inside  the  moment  method 
([z^,]-[a,]  =  [F,]  )  are  calculated  in  the  spectral  domain.  Since  the  problem  is  defined  in  the 

spherical  coordinate  system  we  apply  the  vector-Legendre  transformation  to  the  patch  current  [2] 
-  [4].  This  technique  transforms  the  three-dimensional  problem  into  spectrum  of  one-dimensional 
problems. 

K1 

J(r,0,(|))=  Je  =  J]  (4a) 

_  m=-oo  «=|m| 

j(r,n,m)  = - -  f  fL(«,/n,0)J(r,0,(|»)sin0e'^'^<5?0r/<|)  (4b) 

2nS(n,m)Jl^ 


fjirhm,Q)=  0 
0 


0 

ai;i'^(co60) 

30 

JmPj^^ccsQ) 

sin0 


0 

-jmP^iccsQ) 

sin0 

aj;l'^(cos0) 

30 


S(n,m) 


2n{n  + 1)(«  +  |w|)! 
(2«  +  l)(«-|w|)! 


(4c) 


(4d) 


Figure  4.1.  Patch  antennas  on  spherical  structure,  (a)  rectangular  patch,  (b)  circular  patch. 


Analysis  of  Curved  Microstrip  Antennas _ 15 

Here  /^"'(cosG)  are  the  associated  Legendre  functions  of  the  first  kind.  The  elements  of  the 
moment  method  matrix  Zy,  can  be  expressed  as: 


patch  V 


VhSin0  dQ 


(5) 


where  is  transposed  test  function,  and  G{m,n)  and  j,  (/«,«)  are  the  Green’s  function  and  the 
basis  function  in  spectral  domain.  After  changing  the  order  of  summation  and  integration  we  get 


=-S  Z  S{n,m)rl^^yj{-m,n)G{m,n)3,{m,n) 


(6) 


n>\m\ 


In  the  case  of  the  coax  feeding,  the  radius  of  the  probe  is  usually  a  veiy  small  fraction  of  the 
wavelength  and  the  quantity  rpaich  -  rano  is  usually  small  in  comparison  to  the  wavelength. 
Therefore,  the  probe  is  modeled  as  a  filament  with  a  constant  current 

»•!>)=  ^^2gjjj0  ^  “0/eerf)^  feed)  ^GND  ^  ^  ^  Rpaich,  (^) 

where  and  are  the  0  and  (ji  coordinates  of  the  coax  probe.  The  microstrip  line  can  be 

also  simply  modeled  by  a  filament  with  a  constant  current  placed  close  to  the  edge  of  the  patch 
since  microstrip  line  and  coax  line  have  similar  field  distribution  around  the  feeding  point. 
Furthermore,  if  the  feeding  point  of  the  coax  line  is  close  to  the  patch  edge,  the  input  impedances 
of  the  patch  antennas  fed  by  a  microstrip  line  and  fed  by  a  coax  line  are  very  similar  [5]  -  [7]. 

Using  the  reaction  theorem  we  can  calculate  the  elements  of  the  excitation  vector  for  coaxially- 
fed  patch  antennas  as 


II 

1 

(8a) 

=  Z  Z  J  G(r,m,n)JXm,n)  dr 

(8b) 

m  probe 


4.1.3  Calculation  of  Green’s  Function 

The  Green's  functions  are  determined  in  spectral  domain,  which  was  one  of  the  reasons  for  using 
spectral  approach.  The  GIDMULT  subroutine  calculates  spectral  domain  Green’s  functions  of 
planar,  circular  cylindrical  and  spherical  multilayer  structures.  For  each  input  value  of  the 
spectral  variables  the  subroutine  calculates  the  value  of  the  electric  and  magnetic  field  at  user- 
specified  locations  inside  or  outside  the  multilayer  structure.  The  structure  parameters  are  given 
as  inputs  to  the  subroutine.  The  material  layers  can  have  arbitrary  complex  permittivities  and 
permeabilities. 
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The  solution  procedure  makes  use  of  the  Fourier  transformation  technique.  The  3D  elements  (Fig 
4.2.a)  are  replaced  by  equivalent  or  physical  current  excitations  of  electric  and  magnetic  type, 
which  are  then  Fourier  transformed.  By  the  vector  Legendre  transformation,  the  3D  excitations 
are  transformed  into  harmonic  current  shells  (see  Fig.  4.2.b).  If  the  source  is  infinitely  thin  in  r 
direction,  we  get  one  discrete  current  shell  per  source,  otherwise  we  get  a  continuous  distribution 
of  current  shells  in  r-direction.  The  E-  and  H-fields  induced  by  the  harmonic  current  sources 
have  the  same  harmonic  variations  in  0  and  <!»  as  the  source.  Therefore,  only  the  field  variation 
in  the  direction  perpendicular  to  the  boundaries  is  unknown,  so  we  have  a  harmonic  one¬ 
dimensional  (ID)  field  problem.  In  this  way  the  spectral  domain  problem  is  interpreted  as  a  ID 
spatial  domain  problem  consisting  of  the  ID  multilayer  structure  and  harmonic  ID  sources  in  the 
form  of  current  shells. 

The  harmonic  ID  problem  is  solved  by  making  use  of  equivalent  problems,  one  for  each  layer 
(Figs.  4.2.C  and  4.2.d).  The  unknowns  are  the  tangential  E-  and  H-  fields  at  the  layer  boundaries. 
Since  the  variation  of  the  E-  and  H-fields  in  the  direction  tangential  to  the  boundaries  is 
harmonic  with  known  periodicity,  we  only  need  to  determine  the  complex  field  amplitudes  at  the 
interfaces,  i.e.,  we  have  four  unknowns  per  boundary.  The  algorithm  connects  all  equivalent 
subproblems  into  a  system  of  A{Nlayer-\)  linear  equations  with  the  same  number  of  unknowns 
(Nlayer  denotes  the  number  of  layers).  Once  the  amplitudes  of  the  tangential  fields  have  been 
determined,  it  is  easy  to  determine  the  field  amplitudes  anywhere  in  the  multilayer  structure  by 
applying  the  homogeneous  region  equivalent  of  the  layer  inside  which  we  want  to  determine  the 
field  value.  The  core  problem  in  the  formulation  is  to  calculate  the  E-  and  H-  fields  due  to  a 
harmonic  current  shell  in  a  homogeneous  region  (Fig.  4.2.e).  More  details  about  the  GIDMULT 
algorithm  can  be  found  in  [8]. 


4.1.4  Numerical  Considerations 

4.1.4.1  New  definition  od  vector-Legendre  transformation  and  normalized  Legendre 
polynomials 

We  notice  numerical  problems  in  calculating  Legendre  polynomials  and  S(n,m)  for  large  n, 
specially  if  m  »  n  when  both  Legendre  polynomials  and  S(n,m)  become  very  large,  and  the 
solution  is  unstable  (both  terms  are  needed  for  performing  the  vector-Legendre  transformation; 
notice  that  the  number  of  used  modes  depends  on  the  structure  radius).  To  find  a  numerically 
stable  solution  first  we  have  introduced  an  alternative  definition  of  the  vector-Legendre 
transformation: 

J(r,e,(t))=  ^  -r- ^  L(«,»i,e)  j(f',«,w)  ,  (9a) 

m=-c«  yjz.'K 

n  ji  I  _ 

j(r,«,/n)=  [  f -  -  L(«, m,Q ) J(r,6 ,(|) ) sin9  e~^"^dQ .  (9b) 

i^2nS{n,m) 

The  expressions  for  L(«,w,0)  and  S(n,m)  are  the  same  as  before  (see  equations  (4.c)  and  (4.d)). 
The  new  definition  of  the  vector-Legendre  transformation  and  the  inverse  vector-Legendre 
transformation  is  symmetric,  i.e.  in  both  cases  Legendre  polynomials  and  their  derivatives, 

appearing  in  matrix  L(«,w,0),  are  divided  with  yjs(n,m) . 
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Second  step  in  finding  numerically  stable  solution  has  been  in  defining  normalized  Legendre 
polynomials  and  their  derivatives  as 


pFi  ^  pFI 


l(«-|m)!  dFn  (cos9 )  _  dP}!^ (cos9 )  m )! 

|(«+/k|)!’  59  50  'y  («+/»!)! 


Notice  that  for  normalization  factor  we  have  chosen  the  rapidly  growing  term  inside 
(see  eq.  (4.d)).  By  using  the  following  asymptotic  formulas  [9] 


pW(cos0)  = 


r(«+|w|+l)  2 

r(«  +  3/2)  in  sine 


f  iV  n  |»i|7t 

n  +  -  0 - — — 

[[  2)  4  2 


n->oo  T(n  +  b) 


the  normalized  Legendre  functions  converge  for  large  n  to  the  expression 


limPj'"'(cos0)  =  limPj'”'(cos0) 


i(n+|/«|)!  Vjt«sin0 


(  iV  n  \m\n 

cos  «  +  — 0 — + - 

I  2j  4  2 


Notice  that  it  is  numerically  stable  to  calculate  the  upper  expression.  The  term  L  /  ^S(n,m)  in 
equations  (9a)  and  (9b)  enables  us  to  calculate  normalized  Legendre  polynomials  instead  of 
Legendre  polynomials,  and  2n(n+iy(2n+l)  instead  of  S(n,m),  both  of  them  are  numerically 
stable.  Thus,  divisions  of  very  large  numbers  are  avoided. 

The  recursive  equations  for  the  normalized  Legendre  polynomials  and  their  derivatives  are: 

(z)V(«  +  l+|w|)(«+|w|)  = - r-T77- 

n-\m\+\  (j4) 

•  |(2«  + 1)  z  P^  (z)^(«+|ot|)(«-|ot|+1)  -(n+\m  (z)^in-\m\Xn-\m\+\)  j  , 

— |m| 

(z^  - 1)  ^Jn+  \m\=  nzI^{z)^n-¥  \  m  \  -  («+ 1  m  \)p[-i(z)^Jn-  \m\.  (15) 

5z 

With  the  new  definition  of  vector-Legendre  transformation  the  elements  needed  for  moment 
method  procedure  are  calculated  as: 

^Ji  "  S  n)^{rn,  n)3 ,  {m,  n)  (16) 


m  w=|w|  “y  probe 


(17) 
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4.1.4.2  Asymptotic  formulas  for  spherical  Bessel  and  Hankel  functions 

The  Green’s  functions  for  spherical  structures  contain  the  Schelkunoff  type  of  spherical  Bessel 
and  Hankel  functions,  see  Appendix  II.  We  notice  numerical  problems  in  calculating 
Bessel/Hankel  functions  of  large  order  (e.g.  of  order  larger  than  200;  this  happens  when  the 
radius  of  the  structure  is  large).  Therefore,  special  techniques  are  needed  to  handle  cylinders 
with  large  radii.  We  have  implemented  Debye’s  asymptotic  formulas  for  Bessel  and  Hankel 
functions  of  large  orders  [9],  [10] 


J.ix)  = 


g(m+1/2)(tanh  o  -a )  ^ 


■sj4(m  + 1/2)  tanha 


1  + 


3t  -  5t^ 
24(m  +  l/2) 


a  j 


^(»f+l/2).(a  -  tanh 

^(ot+ 1/2) -tanha 


1  - 


3t  -  5t^ 
lAim  +  yi) 


(18a) 


(18b) 


where  a  and  t  are  defined  as  cosha  =  (w-fl/2)/x  and  t  =  cotha.  Notice  that  the 
exponential  parts  of  Bessel  and  Hankel  approximate  formulas  have  opposite  behavior,  i.e.  the 
arguments  of  the  exponential  functions  have  the  same  absolute  values  and  opposite  signs.  On  the 
other  hand  Green’s  functions  for  spherical  multilayer  structures  can  be  written  in  terms  of 
products  J„(  •  )  •  )  (see  Appendix  II  for  exact  expressions).  Therefore,  Debye’s 

asymptotic  formulas  are  applied  to  products  •/„(•)  ^(  •  )  with  extracted  exponential  parts, 

and  thereby,  we  avoid  numerical  problems.  The  condition  for  order  m  of  Bessel/Hankel 
functions  at  which  asymptotic  formulas  start  to  be  used  take  care  that  l/w-jcl^lxl'^^^for  all  used 
arguments  X  [10]. 


4.1.5  Mutual  Coupling 

First  step  in  determining  iS-parameters  is  to  determine  the  currents  on  all  patches,  i.e.  the  integral 
equation  for  tangential  components  of  the  electric  field  is  solved  by  applying  the  moment 
method.  For  N  patches  in  the  array,  we  need  to  solve  the  linear  system  [2]  [a*]  =  [F*]  A/ times  {k 
=  once  for  each  excitation  port.  In  more  details,  [F*]  vectors  correspond  to  a  physical 

situation  in  which  a  unit  current  is  entering  the  M  port  while  the  remaining  N-1  ports  are  open- 
circuited.  Fortunately,  the  matrix  [Z]  is  unchanged  in  all  cases. 

After  determining  the  amplitudes  of  basis  functions  [a*]  we  calculate  the  voltage  at  port  /  by 

Zl  k 

“/■  >  i.e.,  the  Ik  element  of  the  impedance  port  matrix  is 


(19) 


Scattering  matrix  elements  are  calculated  from  the  impedance  port  matrix: 


(20) 


[SM  )( [Z^^'-'  ]+[z®] )-‘ 

where  [Z^]  is  a  diagonal  matrix  with  elements  2f  -  the  characteristic  impedance  of  the  feeding 
transmission  lines. 

The  elements  of  the  moment  method  matrix  Zji  are  excitation  vector  Vj  are  calculated  in  the 
spectral  domain.  The  expression  for  element  of  the  moment  method  matrix  Zji  is: 

=  "S  Z  jy  n)G(m,  n)j ,  (m,  n)  (21) 

m  w>|m| 

In  matrix  notation  equation  (21)  can  be  written  as: 

[2]^,=-£X[^L’  ZjXn,m)  =  rl^,J^j(-m,n)d(m,n)jXm,n)  (22) 

m  n>\m\ 

In  the  array  case  basis  and  test  function  can  be  located  on  different  patches.  In  that  case 
calculation  of  elements  Z.  is  numerically  complex. 

Let  us  first  consider  mutual  coupling  between  two  patches  displaced  in  the  (|)  direction  (see  Fig. 
4.3).  If  we  assume  that  the  test  function  is  placed  on  the  first  patch  and  the  basis  function  is 
placed  on  the  second  patch,  then  the  element  of  moment  method  matrix  is  given  as  : 

Z)f (n,  m)  =  n)f  G{m, ri)3] (m,  n) 

=  rla,ch  ) 


Figure  4.3.  Patch  array  of  two  patches  displaced  in  the  (|)  direction. 
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Here  the  upper  indexes  1,2  denote  the  patch  number  where  the  basis/test  function  is  located,  and 
the  index  c  denotes  the  centered  basis  or  test  function  (i.e.  the  center  has  coordinates  0  =  nil  and 
(|)  =  0.  In  matrix  notation  the  moment  method  matrix  has  a  shape 


<7CC  y^cc  ycc 

yCCyCC  yCC 


yCCyCC  yCC 
yCCyCC  yCC 

yccycc  ycc 

^2\^22""^2B  gMPj-Pi) 


yccycc  ycc 
_^B\^B2—^BB_ 


ycc  yCC  rycc 
ryccycc  rycc 

^2\^22-—^2B  ^MPi-P2) 


rycc  rycc  rycc 

J^B\^B2""^BB  , 

rycc  rycc  ryCC 

^\\^\2""^\B 

ryCC  rycc  yCC 

^2\  ^22""^2B 


rycc  rycc  ryCC 
_^B\^B2""^BB  _ 


where  B  denotes  the  number  of  basis/test  functions  defined  at  each  patch.  Notice  that 
numerically  there  is  no  difference  in  analyzing  a  single  patch  and  an  array  of  patches  with 
centers  located  at  the  line  0  =  90®. 

For  calculating  mutual  coupling  between  two  patches  spaced  in  0  direction  (see  Fig.  4.4)  it  has 
been  found  convenient  to  introduce  local  coordinate  systems  with  origins  located  at  the  center  of 
each  patch  antenna.  The  vector-Legendre  transformation  of  the  basis/test  function  located  at  the 
displaced  patch  is  calculated  using  transformation  of  coordinates  from  global  to  local  coordinate 
system.  In  this  way  the  basis  function  in  the  spectral  domain  is: 


Xr,n,m)=  y  '  f  L(m,«,0)sin0  e-^'^J(^0 ’,<|) ') 

j2iiS(n,m)  5(0  ,<|)  ) 


where  0  ’  and  (|)  ’  are  coordinates  in  local  coordinate  system  and  |9(0  ,(l))/5(0  ')|  is  the  Jacobian 

determinant 


50 

50 

50' 

5(t>' 

5<|) 

5(j) 

50' 

5(|»' 

The  transformation  from  global  to  local  coordinate  system  is  given  by  [1 1] 


COS0  '  =  sin(a  -Jt  /2)sin0  cos(<t)  -  P)  +  cos(a  -7t  /2)cos0 


_  cos(a  -Jt  /2)sin0  cos((l)  -  P)-sin(a  -Jt  /2)cos0 
sin0  sin((|)  -  P) 
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Here  a  and  p  are  the  0-  and  (j)-coordinates  of  the  patch  center.  To  make  calculations  faster  we 
have  “forced”  separation  of  0  and  (|)  variables  (which  is  strictly  not  possible)  and  approximated 
the  double  integral  with  calculating  numerically  only  0  integration.  In  (t)-direction  we  have 
supposed  that  the  patch  has  constant  width  (i.e  that  the  patch  edges  follow  the  constant  <()-lines), 
and  consequently  we  have  performed  the  (jj-integration  analytically.  Therefore,  for  each  patch 
with  different  0  coordinate  we  need  to  numerically  calculate  the  following  integral 

jE(w,«,0)sin0  J(0 ’)^ci0 ' 

The  moment  method  matrix  for  two  patches  displaced  in  ©-direction  using  B  basis  function  is 
given  as; 


lyCC  7CC 

•yCC 

••^\B 

tyCC  lyCC 

^2X^22" 

7CC 

''^2B 

p^l2^12  pyl2 

yllyXl  yXl 

^21^^22 


rycc  nrcc  nrcc 

n2X^2X  rw2X 
^n^i2—’^iB 
7211721  721 

^21^22  ••••^2B 


7I27I2  7I2 

ryCC^CC  7CC 
Zrjj  Z,j2 

ryCC^CC  7CC 


721 721  721 

^B1^B2****^BB 


rjcc  ycc 
,^B\^B2' 


(28) 


Notice  that  boldfaced  elements  of  the  moment  method  matrix  have  to  be  calculated  using  basis 
function  in  the  local  coordinate  system.  If  the  patches  are  displaced  both  in  0  and  (()  directions, 
then  we  need  to  combine  both  described  methods. 


Figure  4.4.  Patch  array  of  two  patches  displaced  in  the  Q-direction. 
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The  accuracy  of  calculation  of  basis  function  in  spectral  domain  is  shown  in  Fig.  4.5.  The  S21 
coefficient,  calculated  both  in  a  rigorous  and  in  an  approximate  way,  is  compared  for  patches 
displaced  in  E-  and  H-plane.  One  of  the  patches  was  placed  at  center  position  (0  =  tU2,  (|i  =  0), 
and  the  other  patch  was  moved  in  the  0-direction.  The  dimensions  of  the  patches  are  5.1  cm  x 
5.1  cm,  and  the  considered  spherical  structure  has  the  following  parameters:  18.7  cm,  = 

1.0;  h  =  0.52  cm.  It  can  be  seen  that  we  can  apply  the  approximate  way  of  calculating  spectral- 
domain  basis  functions  except  for  rectangular  patches  placed  in  small  regions  around  the  poles. 


Figure  4.5.  Comparison  of  rigorous  and  approximate  way  of  calculating  vector-Legendre 
transformation  of  basis  and  test  functions,  (a)  coupling  in  E-plane,  (b)  coupling 
in  H-plane. 


4.1.6  Far-Field  Calculations 

The  far  field  radiation  pattern  is  obtained  as  follows.  If  we  consider  for  example  the  ^- 
component  of  the  electric  field  in  the  outermost  region  with  the  r-coordinate  larger  than  the  r- 

coordinate  of  the  patch,  we  have  only  outward-traveling  waves  described  by  a^^  H^^ik^r) . 
Therefore,  in  the  outermost  region  we  can  connect  the  (j)-component  of  the  electric  field  with 
different  r-coordinates  ri  and  r2as 


4  ('*2 


(29) 


Here  r/  represents  the  r-component  of  the  far  field  pattern.  The  final  solution  is  obtained  by 
superposing  the  spectral  solutions,  see  eq.  (4).  When  calculating  far  field  of  spherical  array  it  has 
been  found  convenient  to  introduce  a  local  coordinate  system  with  the  origin  located  at  the 
position  of  each  antenna  element.  For  circular  patches  the  z’-axis  of  local  system  is  oriented 
along  the  direction  of  the  pattern  maximum  of  the  individual  element  and  the  x’-axis  is  taken  in 
the  z-z’  plane.  If  the  center  of  the  patch  element  is  at  the  point  P(r,  a„  ,  the  coordinates  in 
the  local  coordinate  system  are  determined  using  the  following  transformation  relations  (see  [11] 
for  details): 
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cos0'=sina„sin0  cos((|)  -  p^)  +  cosa„cos0 

. ,  cosa„  sin0  cos((t)  -  P-„)  -  sina„  cos0 
sin0  sin((|)  -  p„„) 

COS0  sin  a„  cos((|)  -  P„„) - sin0  cosa„  >  ^  sina„sin((|) -p„„) 
sin0 '  ^  sin0 ' 

sina„ sin((|)  -  p„„) ,  cos0  sina„  cos((|)  -  p„„)-sin0  cosa„ 

.  ^  Q  — — -  “ 


sin0 '  °  sin0 ' 

0  and  (|>  field  components  produced  by  antenna  element  placed  at  the  point  P(r,  a„ ,  ^nm)  are: 

>®nPnm 


(30a) 

(30b) 

(30c) 

(30d) 


cos0sina„cos((t>-p„;;,)-sm0cosa„  sma„  sm((t)  -  P«;„)  ^ 

- r-r^ - /0(O  .'P  )  + - r-r, - ><P  ) 

sm0'  sm0' 


""•P  P/i 


(31a) 

sina^sin((|)-p„;;,)^^^^^,  cos0  sina„  cos((|>  -  p„;„)-sin0  cosa^  /y(0',i|i’) 

(31b) 


For  rectangular  patches  the  x  ’  axis  of  local  system  is  oriented  normal  to  the  patch,  i.e.  for  center 
element  axes  x  and  x'  are  in  the  same  direction  (see  Fig.  4.1).  If  a„  and  P„m  are  the  0  and  <|) 
coordinate  of  the  patch  center,  the  coordinates  in  the  local  coordinate  system  are  determined 
using  the  following  equations 

COS0  '  =  sin(a„  -Jt/2)sin0  cos((t)  -  p„„)  +  cos(a„  -n/2)cos0 

^ ^ ,  cos(a„  -7t/2)sin0  cos((|)  -  p„„)-sin(a„  -7i/2)cos0 
sin0sin((|)-p„J 


(32a) 

(32b) 


ee'  =  - 


COS0  sin(a„  -7t/2)cos((|)  -  p„„)-sin0  cos(a„  -n/2)  ^  ,  sin(a„ -7t/2)sin((|)  -  P„„) 


sin0' 


ee  +  - 


sin0' 


(32c) 


e*'  =■ 


sin(a„  -7i/2)sin((|>-p„„)  ^  cos0  sin(a„  -rt/2)cos(<t)  -  P„„,)-sin0  cos(a„  -n/2) 


sin0' 


-ee 


sin0' 


(32d) 


Far  field  produced  by  spherical  antenna  array  is: 

EP.*)' ZEa.|S.„(B'.+') 

n,m 


(33) 
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In  order  to  steer  the  beam  throughout  the  whole  hemisphere  the  arrays  should  have  a  grid 
suitable  for  spherical  structure,  such  as  icosahedron  structure  depicted  in  figure  4.6.  In  this 
example  the  basic  grid  triangle  element  corresponds  to  1/16  of  the  triangle  that  builds  up  the 
icosahedron  structure,  i.  e.  side  of  the  basic  grid  triangle  equals  to  1/4  of  the  side  of  the 
icosahedron  triangle. 


Figure  4.6.  Spherical  array  following  the  icosahedron  grid 


4.1.7  Effects  of  the  Sphere  Curvature 

The  effects  of  the  sphere  size  on  the  input  impedance  and  radiation  pattern  will  be  illustrated  in 
the  following  examples.  The  quasi-rectangular  patch  of  dimension  2.5  cm  x  4.0  cm  is  printed  on 
the  single-layer  dielectric  substrate  with  dielectric  constant  Sr  =  2.52  and  thickness  /i  =  1.58  mm. 
Both  0-  and  (|>-polarized  patches  are  considered,  and  the  patch  length  in  resonance  direction  is 
2.5  cm  in  both  cases.  The  difference  between  two  cases  is  in  the  patch  shape,  i.e.  the  patch  edges 
follow  the  constant  0  and  (|)  lines  (see  Fig.  4.7).  Therefore,  the  0-polarized  patch  has  constant 
length  in  resonance  direction  (2.5  cm),  while  the  resonant  length  of  the  (|)-polarized  patches 
varies  from  2.5  cm  in  the  center  of  the  patch  up  to  1.74  cm  at  the  patch  edge  (for  2.5  cm  sphere 
radius).  The  feed  point  is  placed  1.1  cm  from  the  patch  center. 


Figure  4. 7.  Geometry  of  the  problem  -  azimuthal  equatorial  projection. 
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The  influence  of  the  sphere  radius  on  input  impedance  of  the  spherical  patch  is  shown  in  Figures 
4.8  and  4.9.  For  comparison,  the  measured  input  impedance  of  the  planar  patch  antenna  with 
same  dimensions  is  also  given.  It  can  be  seen  that  the  sphere  radius  mostly  influences  the 
resonant  frequency,  and  the  change  of  resonant  frequency  is  stronger  for  the  (|)-polarized  patch 
(since  the  average  patch  length  is  shorter  for  the  {t)-polarized  patch).  The  resonant  resistance  is 
almost  constant  for  the  0-polarized  patch  and  it  is  slowly  changing  for  the  (t)-polarized  patch. 

The  effect  of  the  radius  of  the  ground  plane  (grounded  shell)  on  radiation  pattern  is  given  in  Fig. 
4.10.  For  comparison,  the  calculated  radiation  pattern  of  the  planar  patch  antenna  with  same 
dimensions  is  also  given.  The  operating  frequency  is  3.6  GHz.  It  can  be  seen  that  with  enlarging 
the  radius  of  the  sphere  the  main  lobe  of  the  spherical  patch  antenna  approaches  the  main  lobe  of 
the  planar  counterpart.  As  expected,  the  back-radiation  is  smaller  for  spheres  with  larger  radius. 

For  building  up  a  spherical  array  grid,  the  icosahedron  structure  is  very  suitable  basis.  Beside  the 
beam  steering  throughout  the  whole  hemisphere,  this  grid  enables  beam  width  variation  by 
changing  the  number  of  excited  elements.  To  illustrate  this,  radiation  pattern  calculation  of  a 
spherical  array  with  a  grid  as  depicted  in  figure  4.6  (where  basic  triangle  side  is  1/4  of  the 
icosahedron  triangle)  is  performed.  Activated  area,  i.  e.  amount  of  excited  elements  was  varied 
from  angle  a  =  15,  30  and  45  degrees.  This  corresponds  to  6,  16  and  31  activated  elements, 
respectively.  The  example  in  Fig.  4.6b  corresponds  to  a  =  30  degrees.  The  array  grounded  shell 
radius  is  20  cm,  size  of  the  patches  is  2.5  cm  x  2.5  cm,  and  substrate  and  frequency  of  analysis 
are  the  same  as  in  Figs.  4.8-4.10.  The  patches  are  phase  corrected  in  order  to  compensate  the 
propagation  path  difference.  The  radiation  pattern  for  all  three  different  activated  areas  is  given 
in  figure  4.11.  As  expected,  with  enlarging  the  number  of  activated  elements  the  main  beam  is 
narrower.  However,  the  back  radiation  increases  with  increase  of  number  of  activated  elements. 

The  influence  of  the  sphere  radius  on  mutual  coupling  level  is  shown  in  Fig.  4.12.  Two 
rectangular  patches  of  dimension  10.57  cm  x  6.55  cm  are  printed  on  the  single-layer  dielectric 
substrate  with  dielectric  constant  =  2.5  and  thickness  h  =  \  .575  mm.  The  position  of  the  probe 
is  1.5  cm  from  the  center  of  the  patch.  The  working  frequency  for  planar  patch  is  1.41  GHz.  For 
spherical  case  the  working  frequency  is  changed  a  bit  in  order  to  obtain  good  matching  (the 
needed  shift  of  resonant  frequency  depends  on  the  sphere  radius).  The  magnitude  of  the  5'2i 
parameter,  both  in  E-  and  H-plane,  is  shown  as  a  function  of  spacing  between  patch  edges.  For 
comparison,  the  measured  results  of  planar  case  are  also  given  [12].  It  can  be  seen  that  the 
coupling  is  weaker  for  spherical  case  in  comparison  to  the  planar  one,  especially  in  H-plane. 
However,  for  smaller  structure  radii  mutual  coupling  starts  to  be  oscillatory. 


normalized  impedance  -  imaginary  part  normalized  impedance  -  real  part 
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frequency  (GHz) 


(b) 


Figure  4.8.  Input  impedance  of  a  Q-polarized  quasi-rectangular  patch  printed  on  a  spherical 
structure  for  different  sphere  radii;  (a)  real  part,  (b)  imaginary  part. 
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(a) 


frequency  (GHz) 
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Figure  4.9.  Input  impedance  of  a  ^-polarized  quasi-rectangular  patch  printed  on  a  spherical 
structure  for  different  sphere  radii;  (a)  real  part,  (b)  imaginary  part. 
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(b) 


Figure  4.10.  Radiation  pattern  of  a  spherical  rectangular  patch  antenna  for  different  sphere 
radii;  (a)  H-plane,  (b)  E-plane. 
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(a) 


(b) 

Figure  4.11.  Radiation  pattern  of  a  spherical  patch  array  with  elements  placed  uniformly  by 
utilizing  the  symmetry  properties  of  the  icosahedron,  (a)  H-plane,  (b)  E-plane. 
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Figure  4. 12.  Magnitude  of  S21  parameter  of  a  two-patch  array  as  a  function  of  spacing 
between  patch  edges,  (a)  E-plane,  (b)  H-plane. 
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4.1.8  Laboratory  Model 

In  order  to  validate  the  calculation  results  a  laboratory  model  shown  in  figure  4.13  was 
developed,  which  enabled  simple  measurements  of  different  array  configurations.  The  laboratory 
model  was  built  from  a  copper  half-sphere  of  radius  a  =  18.7  cm  at  which  patches  can  be 
mounted  at  arbitrary  positions.  Mounting  of  a  patch  at  desired  position  is  very  simple:  one  only 
needs  to  place  a  SMA  connector  at  selected  position  on  the  grounded  shell. 

coaxially  fed 
patch 


styrofoam  spacers 
copper  half-sphere 


Figure  4.13.  Sketch  and  photo  of  the  developed  laboratory  model. 

The  selected  dimensions  of  the  patches  are  5.1  cm  x  5.1  cm,  and  the  appropriate  styrofoam  cubes 
are  used  as  spacers  between  patches  and  grounded  shell  (therefore  Sr  =  1 .0;  the  cube  thickness  of 
0.52  cm  is  selected).  Fig.  4.14  shows  the  comparison  of  calculated  and  measured  input 
impedance.  The  mutual  coupling  of  two  patches  in  E-plane  is  studied  as  a  function  of  the 
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frequency  and  spacing  between  the  elements  (Figure  4.15).  The  calculated  and  measured  mutual 
impedances  are  in  a  good  agreement. 


Figure  4.14.  The  comparison  of  calculated  and  measured  input  impedance;  (a)  real  part,  (b) 
imaginary  part. 
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Figure  4. 15.  The  comparison  of  calculated  and  measured  mutual  coupling  coefficient  of 
patches  in  E-plane  as  a  function  of frequency  and  spacing  between  elements; 

(a)  10  cm ,  (b)  20  cm  . 
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As  basis  functions  we  have  used  entire-domain  basis  functions,  i.e.,  basis  functions  that  are 
defined  on  the  whole  patch.  For  rectangular  patches  basis  functions  are  sinusoidal  in  the  current 
direction  and  they  are  of  constant  distribution  in  the  perpendicular  direction,  which  is  a  common 
approximation  since  patch  antennas  are  resonant  structures.  For  example,  the  kth  mode  of  0- 
polarized  current  is 


^^patch^  p 


-Sin 


kn 


20. 


0+0. 


(34) 


and  its  vector-Legendre  transformation  is: 


-  sin(w(t)  ) 
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The  kth  mode  of  (|)  polarized  rectangular  patch  is: 
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•^♦(0»<1>)  =  T - ^sin 


^^paich^p 


and  its  vector-Legendre  transformation  is 
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If  circular  patch  is  analyzed,  then  the  basic  current  mode  is  given  as: 


.  ^  1 

J  =a„7,(a„r^,,,0)cos(|)0 - -  J,(anV.Ae)sin^ 

^polclP 


(38) 


where  ay;  denotes  constant  (an  rpaich%  =  1.84118),  and  Jj  is  the  Bessel  function.  The  vector- 
Legendre  transformation  is 


(n,  m)  ■ 


2n  S(n,m) 


\Pn  (COS0)  t/G  +  J  (Xu  sin 0  )  P 
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APPENDIX  II: 

Electromagnetic  Field  Radiated  by  an  Current  Shell 

The  core  problem  inside  the  GIDMULT  algorithm  is  to  calculate  the  electromagnetic  field  in  the 
homogeneous  space  radiated  by  an  electric  current  shell  of  radius 


a)  from  J^in,m) 

-  ,  ,  k  nj.  ^  jHfikrjldcr) 

E,(,r,n,m)  = - -J^in,m)  •  \ 


E,(r,n,m)  = - f^nin  +  l)Jg(n,m)  ■  \  -  . 

coe  r  I 
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b)  from  J.(n,m) 
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HXr,n,m)  =  + 


[  l(krJHi^\kr) 


(41c) 


Here  k  is  the  wave  number.  The  fields  excited  by  the  magnetic  shell  can  be  easily  determined  by 
applying  the  duality  concept.  The  fields  caused  by  the  r-directed  sources  are  evaluated  by  using 
the  transverse  replacement  currents  (see  [4]  for  details). 
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4.2  Analysis  of  aperture  coupled  patch  antennas 
on  circular-cylindrical  and  spherical  structures 


4.2.1  Introduction 

Microstrip  antennas  can  be  simply  fed  by  microstrip  or  coaxial  transmission  lines.  The  microstip 
line  feeding  has  advantage  of  easy  fabrication,  but  it  may  not  be  desirable  due  to  the  radiation  of 
the  microstrip  line.  Coaxial  line  is  not  attractive  from  the  fabrication  point  of  view.  Both  types  of 
feeding  structures  lead  to  narrow-band  patch  designs  (typically  1-3%).  The  aperture-coupled 
feeding  structure  has  advantage  of  both  fabrication  simplicity  and  broadband  patch  antenna 
design  possibility.  Furthermore,  it  has  potential  in  millimeter-wave  applications  and  in 
monolithic  microwave  integrated  circuits  applications. 

The  purpose  of  this  chapter  is  to  present  an  electromagnetic  model  that  enables  the  analysis  of 
aperture-coupled  patch  antennas  mounted  on  a  circular-cylindrical  or  spherical  structure.  The 
analysis  method  is  based  on  the  equivalence  principle  which  is  used  to  separate  the  inner 
problem  (feeding  structure  inside  the  grounded  tube  or  grounded  shell)  and  the  outer  problem 
(patch  antenna  itself).  The  reciprocity  theorem  is  used  to  find  the  interaction  between  the 
feedline  and  the  aperture.  The  antenna  structure  is  rigorously  taken  into  account  by  using  proper 
Green’s  functions,  and  the  array  is  fully  analyzed  by  means  of  element-by-element  approach. 
The  integral  equation  is  numerically  solved  using  spectral  domain  moment  method.  The 
developed  program  calculates  radiation  pattern  of  the  antenna,  input  impedance  at  each  port  and 
mutual  coupling  between  the  elements  of  the  array. 

4.2.2  Method  of  Analysis 

The  analysis  of  microstrip  patch  arrays  on  circular  cylindrical  structures  is  given  in  [13].  Here 
we  will  give  only  the  analysis  model  of  aperture-coupling,  i.e.  we  will  extend  the  analysis  to  the 
new  feeding  structure.  Furthermore,  the  same  model  can  be  applied  both  to  circular-cylindrical 
and  spherical  structures.  Therefore,  we  will  give  the  description  of  cylindrical  case  only  since  the 
extension  of  formulation  to  spherical  problems  is  straightforward. 

4.2.2.1  Analysis  of  one  aperture-coupled  patch 

First  we  will  describe  the  feeding  model  for  single  aperture-coupled  patch  antenna.  The 
geometry  of  the  patch  is  given  in  Figs.  4.16  and  4.17  (notice  that  the  Figure  4.17  shows  only  a 
small  part  of  a  cylindrical  or  spherical  antenna).  The  analysis  of  aperture  feeding  based  on 
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reciprocity  theorem  was  developed  for  planar  patches  by  D.  Pozar  [14].  Here  we  will  give  the 
formulation  for  patches  on  circular-cylindrical  structures. 


It  is  assumed  that  a  quasi-TEM  mode  is  propagating  along  the  microstrip  line  oriented  in  z- 
direction 

etp.d))  =  e„p  -l-exA 


Figure  4.16.  Aperture-coupled  patch  antenna,  (a)  Circular-cylindrical  patch  antenna,  (b) 
Spherical  patch  antenna. 


of  length /s/wb 


Figure  4. 1 7.  Aperture-coupled  patch  antenna.  The  antenna  is  mounted  on  cylindrical  or 
spherical  structure. 
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(43) 


Here  p  is  a  propagation  constant,  and  e  and  h  are  determined  from  adequate  Green's  functions. 
The  aperture  represents  a  discontinuity,  therefore  we  can  write 


for  z  <0 

(44a) 

]  ’ 

]  _ 

[TE^  , 

for  z  >  0 

\h*+rh-, 

for  z  <  0 

(44b) 

[TH^  , 

for  z>0 

where  R  and  T  represent  reflection  and  transmission  coefficients.  By  applying  the  reciprocity 
theorem  [14]  on  the  total  fields  E,  H  and  on  the  forward  propagating  fields  E*,H*  we  get: 

jExH^  ■ds  =  jE^  xH  ds,  (45) 

s  s 

where  S  represents  the  closed  area  consisting  of  three  parts: 

-  S(,=  effective  cross-section  of  microstrip  line  (0  <  p  <  Pgnd.  0^  ^<2n) 

-  S„=  aperture  area 

-  S„  =  walls  of  microsrip  line  at  the  PEC  grounded  tube  (p  =  pgnd,  0^  ^<2%) 


On  we  have  hxE  =nxE*  =0,  thus  there  is  no  contribution  when  integrating  over  .  On 
fixE*  =0,  but  =  Fo«x(f  e^,  where  is  the  unknown  field  distribution  in  the 
aperture,  and  we  will  represent  it  via  piecewise  sinusoidal  (PWS)  basis  function: 


V  e" 


■K 


smk^(L/2-  I 


'GND 


W sin  EL/ 2 


for  r  <SW/2, 


'CND 


<SL/2. 


(46) 


Here  Ice  is  the  effective  wave  number  of  PWS  basis  function,  and  is  the  unknown  amplitude  to 
be  determined.  For  Ice  it  is  convenient  to  use  the  mean  value  of  wave  numbers  of  two  surrounding 
layers  (the  layers  next  to  the  grounded  PEC  tube  that  are  inside  and  outside  the  tube).  Equation 
(32)  can  be  written  in  the  following  form: 

^ (x,  y)h  (^’ 

^  s,  ^ 

For  narrow  slots  there  is  no  need  to  include  term  in  equation  (34)  since  its  influence  is 
small.  By  applying  the  reciprocity  theorem  once  more 
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\ExH-  ■ds  =  \E-xH-c^,  (48) 

s  s 

we  get  expression  for  transmission  coefficient  T : 

r  =  1  -  \e:ix,y)h^ {x,y) ds  =  \-R.  (49) 

Since  we  have  three  unknowns  R,  T  and  we  need  one  more  equation.  We  will  apply  the  PEC 
equivalence  theorem  [14]  and  by  this  we  will  separate  the  inner  and  the  outer  region.  Therefore, 
we  will  introduce  the  equivalent  magnetic  current  at  the  aperture  =-«xE.  The  third 
equation  is  a  consequence  of  the  continuity  of  the  magnetic  field  across  the  aperture: 

h;=hI+h[,  (50) 


where 

=  outside  magnetic  field  (p  >pc/®)  due  to  the  aperture  field  V^el 
H'^  =  inside  magnetic  field  (p  <Pgm))  due  to  the  aperture  field 
H(  =  inside  magnetic  field  (p  <poM3)  due  to  microstrip  line  quasi-TEM  mode 

At  z  =  0"  (or  z  =  0^  since  \  -R  =  T),  it  is  valid: 

Hl={\-R)h^.  (51) 


Therefore, 


Hl-H[  =V,  |(G;f"((t),z;(|,o,z,)-G;f"((|),z;<t>„,Zo)).<(<|)o,2o)^f^,.  (52) 


Here  and  are  respectively  the  Green’s  functions  for  external  and  internal  regions, 

A 

and  <|)  •  e"  ((]) ,  z)  is  the  equivalent  magnetic  current.  Combining  equations  (50)  -  (52)  we  get 


K  j(G;f"((|),z;(j)o,Zo)- 


G'I,hm 
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After  multiplying  (40)  with  e"  and  integrating  over  we  get 


SaS, 


=  (l-R)le:(<^,z)h,(<^,z)ds 


(54) 


or 


Analysis  of  Curved  Microstrip  Antennas _  42 

Fo(r  +  7')  =  Av(l-i?),  (55) 

r  =  j  ■  <((t)o,Zo)  ds  ds,  (56) 

SaSa 

r' =  JJ<((|),z)G;f*"((l),z;(t>o,2o)-<(<l>o.^o)^'S^So  (57) 


When  the  patch  is  present  we  need  to  add  contribution  of  the  patch  as  well  (in  order  to  fulfill  the 
continuity  of  magnetic  field).  For  that  first  we  need  to  determine  the  current  distribution  on  the 
patch,  which  is  done  by  applying  the  moment  method.  The  patch  current  is  expanded  in  entire 
domain  basis  functions  multiplied  by  unknown  coefficients  a/ 

J((t),z)  =  |;a,J,((l),z).  (58) 

/=! 

We  formulate  the  electric  field  integral  equation  (EFIE)  by  enforcing  the  boundary  condition  that 
the  total  tangential  electric  field  must  be  zero  on  the  surface  of  the  patch.  The  equation  is 
satisfied  in  an  average  sense  by  weighting  the  fields  with  test  functions  (that  are  equal  to  basis 
functions  -  Galerkin  method).  We  get  the  matrix  equation  to  solve 

[z]kHy],  (59) 

where  [a]  is  vector  of  unknown  coefficients  a/,  and  impedance  matrix  [z]  and  voltage  vector 
\v]  are  calculated  in  the  following  way: 

Z„„  =  -  J  Jj,(<|),z)  •  ;(|)o,Zo)  •  Ji((|),z)^fs  dso,  (60a) 

K  =  j  jj/<|>,z)-G^"((|),z  ;(|)o,Zo)  4 <(<l>o.^o)  ^  dso-  (60b) 

The  weighted  magnetic  field  at  the  aperture  due  to  the  patch  current  is  given  by  (notice  that  we 
have  used  the  reaction  theorem  in  the  derivation): 

Y"  =  [V]'[ZY[V]  (61) 

Therefore,  the  modified  equation  (55)  with  added  weighted  magnetic  field  from  the  patch  is 

Fo(y“  +Y‘  +  r)  =  Av(l  -  R).  (62) 

Combining  equations  (62)  and  (47)  we  get 

Av^ 

_ 

Av"  +  2(r +r +r) 


(63) 
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Equation  (49),  T  =  \-R,  actually  means  that  the  slot  represents  a  series  impedance  that  is 
connected  to  the  open  stub  of  length  Istub  (see  Fig.  4.18).  Finally,  we  get  expression  for  input 
impedance  of  one  aperture-coupled  patch: 

^  ^  °  r  r  ~ 

where  Zq  is  the  characteristic  impedance  of  microstrip  line. 


o- 


Zo 


oo 


Zo 


o 


o - oo - o 

Aperture  equivalent  Open  stub  for 

circuit  impedance  matching 


Figure  4. 18.  Equivalent  circuit  of  aperture  in  the  ground  plane 


4.2.2.2  Analysis  of  the  array  of  aperture-coupled  patches 

The  array  of  aperture-coupled  patches  is  analyzed  using  the  equivalent  circuit  of  an  aperture- 
coupled  patch  array  considered  as  a  multiport  device.  Notice  that  each  patch  is  considered  as  two 
port  device  and  that  one  port  is  connected  to  the  open  stub  of  length  ktub  (see  Figs.  4.17  and 
4.18).  If  we  would  consider  one  aperture-coupled  patch  in  this  way,  then  its  S-parameters  will  be 
Sn=S22=7?,  Si2=S2i=7^7-i?.  Using  standard  algebra  with  S-parameters  we  get  the  reflection 
coefficient  Fin  of  the  two  port  device  with  a  load  connected  to  the  second  port  (the  connected 
load  is  characterized  via  reflection  coefficient  Fstub) 


^/n  ~  *^1 1 


1  — 


which  is  the  same  expression  as  (64). 


(65) 


If  we  organize  the  equivalent  2N-port  device  (N  is  number  of  patches  in  the  array)  in  the  way 
that  first  N  ports  are  ports  connected  to  generators  and  the  second  N  ports  are  loaded  with  open 
stub  of  length  then  we  can  divide  the  S  matrix  of  the  2N-port  device  as 


[S.]\ 


(66) 


where  [U ]  is  the  unit  matrix.  Similarly  to  the  one  patch  case  we  get  ,  which  is  equal  to  the 
S-matrix  of  loaded  2N-port  device  with  N  open  stubs,  as 


I 
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[■S'2.]-  (67) 

<  ^  stub  J 

The  elements  of  the  [i",,]  matrix  are  determined  in  the  following  way.  For  each  slot  we  can  write 
(2N-port  equivalent  of  equations  (47)  and  (62)) 

+n  +  yl)  =  Av(5,  -  5,),  (68) 

where  the  index  ij  means  that  the  yth  port  is  excited  and  we  look  at  the  S-parameter/amplitude  of 
the  aperture  field  at  ith  patch,  and  8,y  is  the  Kronecker  symbol.  Once  we  know  the  S-matrix  of  the 
loaded  2N-port  device,  we  determine  all  needed  parameters  in  the  standard  way.  For  example, 
the  impedance  port  matrix  [z'®"]  is  determined  as: 

[ZPor,^  =  [Z‘’]([[/]-[^]  J([t/]  +  [5]^)-'  (69) 

where  [Z^]  is  a  diagonal  matrix  with  elements  Z^,  the  characteristic  impedance  of  the  feeding 
transmission  lines. 
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5  USER'S  MANUAL 


Analysis  of  Curved  Microstrip  Antennas 


46 


5.1  Introduction  to  SMiPA 


5.1.1  Introduction 

This  chapter  gives  the  user’s  manual  for  the  program  SMiPA.  The  abbreviation  SMiPA  means 
the  program  for  analyzing  Spherical  Microstrip  Patch  Arrays.  Since  there  are  different 
requirements  on  calculation  of  radiation  pattern  and  input  impedance/mutual  coupling,  we  have 
actually  developed  two  programs  SMiPA  and  SMiPApat.  The  main  reason  for  splitting  the 
problem  into  two  parts  is  that  radiation  pattern  is  usually  calculated  only  for  one  frequency, 
while  input  impedance  and  mutual  coupling  is  usually  calculated  for  frequency  band  of  interest 
(and  thus  even  hundreds  of  frequency  points  can  be  calculated).  Therefore,  program  SMiPA 
calculates 


•  current  distribution  at  each  patch  in  the  array 

•  input  impedance  at  each  input  port  in  the  array 

•  mutual  coupling  between  each  two  patches  in  the  array 


and  program  SMiPApat  calculates 

•  radiation  pattern  of  the  array  without  taking  mutual  coupling  into  account 

•  radiation  pattern  of  the  array  when  all  mutual  couplings  are  taken  into  account. 


Program  is  written  in  FORTRAN  90,  and  therefore  can  be  compiled  and  run  at  each  machine 
where  FORTRAN  90  is  installed  or  which  support  the  executable  version  of  the  program.  The 
communication  with  the  program  is  made  via  input/output  ASCII  files.  For  case  of  simplicity, 
both  programs  SMiPA  and  SMiPApat  are  using  the  same  input  file  (smipa.in).  The  detailed 
description  of  input  and  output  files  are  given  in  this  chapter. 

We  have  also  developed  a  graphical  user  interface,  called  SMiPAwin,  for  IBM  PC  compatible 
computers  under  Windows  operating  system  (the  description  of  SMiPAwin  is  given  in 
Appendix).  By  this  it  is  easier  to  define  the  input  parameters,  and  to  obtain  graphical 
presentation  of  the  results.  Programs  SMiPA  and  SMiPApat  are  working  independently  of 
graphical  user  interface,  i.e.,  graphical  user  interface  is  used  only  for  writing  the  input  file, 
running  the  SMiPA  and  SMiPApat  programs,  reading  the  output  files,  and  displaying  the  results. 
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5.1.2  Problem  Domain  Description 

The  SMiPA  and  SMiPApat  programs  analyze  microstrip  patch  antennas  or  antenna  arrays  placed 
on  spherical  structures,  as  illustrated  in  Fig.  5.1.  If  an  array  is  analyzed,  following  assumptions 
are  made: 

•  all  antenna  elements  are  identical,  both  in  terms  of  dimensions,  shape  and  feed/port  types 

•  the  patch  elements  are  of  circular  or  rectangular  shape 

•  the  dielectric  structure  can  be  multilayered  with  losses 

•  the  patch  array  must  be  placed  at  the  top  of  one  such  layer 

•  three  types  of  feeding  structures  can  be  analyzed:  coax  feeding,  microstrip  line  feeding  and 
aperture-coupled  patch  elements 

•  multiple  ports  can  be  defined  on  patch,  e.g.  when  designing  circularly  polarized  arrays. 


By  means  of  SMiPA  program  current  distribution  at  each  patch  in  the  array  can  be  calculated,  as 
well  as  input  impedance  at  each  input  port  and  mutual  coupling  between  each  two  patches.  The 
SMiPApat  program  calculates  then  radiation  pattern  of  the  array,  with  or  without  taking  mutual 
coupling  into  account. 


Figure  5. 1.  The  geometry  of  the  spherical  patch  array 
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5.1.3  Text  File  Interface 


The  programs  SMiPA  and  SMiPApat  interact  with  the  user  and  with  each  other  using  ASCII  text 
files.  This  scheme,  with  all  files  in  question  is  presented  in  Fig.  5.2. 


SMiPA. out 


OUTPUT  DATA  | 

*  sij.dat 


Figure  5.2.  Input  and  output  files  for  the  programs  SMiPA  and  SMiPApat. 


The  program  can  calculate  the  following  antenna  characteristics: 

(a)  current  distribution  at  each  patch  in  the  array 

(b)  input  impedance  of  each  patch  in  the  array 

(c)  mutual  coupling  between  each  two  patches  in  the  array 

(d)  radiation  pattern  of  the  array  when  all  mutual  couplings  are  taken  into  account 

(e)  radiation  pattern  of  the  array  without  taking  mutual  coupling  into  account  (fast  calculations 
of  the  radiation  pattern  needed  for  making  the  first  design  of  the  array). 
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The  name  of  the  input  file  for  both  SMiPA  and  SMiPApat  programs  is  SMiPAJn.  Each  line 
possesses  the  value  of  several  variables  and  short  description  of  these  variables.  The  input 
variables,  which  must  have  values  when  the  call  of  the  SMiPA  and  SMiPApat  is  made,  are: 

First  line: 

Fmin.  Fmax.  Nfreauencv.Modecal _ REAL.  REAL.  INTEGER.  CHARACTER 

Fmin  and  Fmax  are  the  minimum  and  maximum  frequency  (in 
GHz).  Nfrequency  is  the  number  of  frequency  points  for  which  the 
array  is  analyzed.  These  values  are  used  only  when  the  call  of  the 
SMiPA  program  is  made.  SMiPApat  program  calculates  the 
radiation  pattern  only  for  one  frequency,  and  that  frequency  is 
Fmin.  If  we  first  calculate  the  patch  currents  at  Nfrequency 
frequency  points  (by  calling  SMiPA  program),  then  we  can  choose 
the  frequency  point  for  which  we  want  to  calculate  the  radiation 
pattern  by  choosing  Fmin.  The  only  condition  is  that  the  patch 
currents  are  calculated  for  the  chosen  frequency  point.  Modecal 
represents  the  calculation  mode  -  R  means  accurate  (rigorous) 
mode  suggested  for  final  design,  and  F  means  less-accurate  (fast- 
design)  mode  suggested  for  first  design. 

Second  line: 

Rgnd _ REAL 

Rgnd  is  the  radius  of  the  grounded  tube  (in  cm). 


Third  line: 

Nlaver _ INTEGER 

Nlayer  is  the  number  of  dielectric  layers. 


Next  Mayer  lines: 

h.  EpsilonR.  Tdelta _ REAL 

Each  line  contains  three  parameters:  h,  EpsilonR  and  Tdelta  which 
describe  the  thickness,  the  relative  permittivity  and  loss  tangent  of 
each  dielectric  layer,  respectively. 


Next  line: 

Noatl .  NDat2.  Lattice  INTEGER.  CHARACTER 

Lattice  represents  the  type  of  the  grid  lattice  -  R  means  rectangular 
lattice,  /  means  icosahedron  lattice,  and  F  means  that  the  position 
of  each  patch  is  read  at  the  end  of  the  input  file.  For  rectangular 
lattice  Npatl  and  Npat2  are  the  number  of  patches  in  the  0-  and  (|)- 
directions,  respectively.  For  icosahedron  lattice  Npatl  is  the 
number  of  rings  in  the  array,  and  Npat2  denotes  the  ratio  between 
the  sides  of  equilateral  triangles  of  the  icosahedron  structure  and  of 
the  considered  array  (for  example,  in  Fig.  4.6  Npatl=3  and 
Npat2=4).  If  the  patch  position  is  read  from  the  file,  then  Npatl 
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and  Npat2  have  no  meaning  and  therefore  they  can  have  arbitrary 
values. 


Next  line: 

Noatlaver _ INTEGER 

Npatlayer  is  the  number  of  the  dielectric  layer  at  top  of  which  the 
patches  are  placed.  In  other  words,  patches  are  placed  between 
dielectric  layers  Npatlayer  and  Npatlayer+\.  If  the  patches  are 
placed  in  the  middle  of  some  dielectric  layer,  an  additional 
dielectric  layer  should  be  introduced,  i.e.  the  considered  dielectric 
layer  should  be  split  into  two  parts  and  the  patches  should  be 
placed  at  the  interface  between  these  two  dielectric  layers. 


Next  line: 

Wpatth.  Wpatphi.  Patshape  REAL.  CHARACTER 

For  rectangular  patches  Wpatth  and  Wpatphi  denote  the  patch 
width  (in  cm)  in  the  0-  and  (j)-  directions,  respectively.  For  circular 
patches  Wpatth  denotes  the  patch  diameter  (in  cm).  Patshape 
represents  the  type  of  the  patch  shape:  R  means  rectangular  shape 
and  C  means  circular  shape. 

Next  line: 

Modefeed _ CHARACTER 

Modefeed  denotes  the  type  of  the  patch  feeding.  M  or  m  means 
microstrip  line  feeding,  C  or  c  means  coax  feeding,  and  J  or  a 
means  aperture  coupled  patch  elements. 


Next  line: 

Nport _ INTEGER 

Nport  denotes  the  number  of  ports  per  patch  (i.e.  how  many  ports 
has  each  patch). 


Next  Nport  lines: 

Feedth.  Feedphi _ REAL 

Feedth  and  Feedphi  determine  the  position  of  the  feeding  line  (in 
cm)  in  the  0-  and  (|)-direction,  respectively.  For  aperture  coupled 
patch  elements  Feedth  and  Feedphi  denote  position  of  the  center  of 
the  feeding  aperture.  Each  line  describes  one  port  of  the  patch. 
This  position  is  relative  to  the  patch  center,  i.e.,  the  feeding  line 
placed  at  patch  center  has  coordinates  Feedth  =  0  and  Feedphi  =  0. 
If  the  patch  is  fed  by  the  microstrip  line,  one  value  should  be 
±Wpatthl2  or  ±Wpatphi/2,  depending  at  which  edge  the  patches  are 
excited.  Feedth  and  Feedphi  also  determine  the  polarization  of  the 
patch  elements  when  SMiPApat  is  used  in  the  simple  mode.  In 
more  details,  when  Feedth  =  0  current  on  patches  is  (|)-directed, 
and  when  Feedphi  =  0  current  on  patches  is  0-directed.  For 
aperture-coupled  patches  the  dimensions  of  the  slot  are  also  used 
in  determining  the  polarization  of  patch  elements.  Otherwise,  the 
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current  on  patches  is  a  linear  combination  of  0-  and  (l»-directed 
basic  basis  functions. 


Next  line: 

Modeexci _ CHARACTER 

Modeexci  determines  if  the  ports  have  equal  excitations,  both  in 
amplitude  and  phase.  F  or  y  means  that  the  ports  are  equally 
excited,  Norn  means  that  the  excitations  are  read  at  the  end  of  the 
input  file  SMiPA.in. 


Next  line: 

Dpatth,  Dnatphi _ REAL 

In  the  case  of  rectangular  array  grid  (i.e.  if  the  variable  Lattice  is 
equal  R)  Dpatth  and  Dpatphi  are  the  distances  between  centers  of 
the  patches  (in  cm)  in  the  0-  and  ([(-directions,  respectively.  If  the 
variable  Lattice  is  equal  /  or  F,  Dpatth  and  Dpatphi  can  be 
arbitrary  numbers 

Next  tine  for  coax  feeding  or  microstrip  line  feeding: 

ZO _ REAL 

ZO  is  the  characteristic  impedance  of  the  input  transmission 
line  (either  microstrip  or  coax  line). 

Next  several  lines  for  aperture-coupled  patch  elements: 

Next  line: 

Slotth.  Slotphi _ REAL 

Slotphi  and  Slotz  are  the  widths  of  the  slots  (in  cm)  in  the 
<[)-  and  z-directions,  respectively. 

Next  line: 

Nlaverin _ INTEGER 

Nlayerin  is  the  number  of  dielectric  layers  inside  the  sphere. 


Next  Nlayer  lines: 

h.  EpsilonR.  Tdelta _ REAL 

Each  line  contains  three  parameters:  h,  EpsilonR  and  Tdelta 
which  describe  the  thickness,  the  relative  permittivity  and 
loss  tangent  of  each  dielectric  layer  inside  the  sphere, 
respectively.  First  considered  layer  is  the  layer  next  to  the 
grounded  shell  that  contains  the  aperture. 


Next  line: 

Nfeedlaver _ INTEGER 

Nfeedlayer  is  the  number  of  the  dielectric  layer  (starting 
from  grounded  shell)  at  top  of  which  a  microstrip  feed  line 
is  placed.  In  other  words,  the  microstrip  feed  line  is  placed 
between  dielectric  layers  Nfeedlayer  and  Nfeedlayer+\ . 
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Next  line: 

Wfeed.  Lstub _ REAL 

Wfeed  is  the  width  of  the  input  microstrip  line  (in  cm). 
Lstub  is  the  length  of  the  open  stub  (in  cm). 

The  next  six  lines  are  needed  only  when  the  call  of  the  SMiPApat  program  is  made.  They  are: 


Phimin.  Phimax.  Nphi _ REAL.  REAL.  INTEGER 

The  minimum  and  maximum  angles  of  (j)  (in  deg)  in  the  azimuthal 
plane,  and  the  number  of  points  which  are  calculated  in  the 
azimuthal  pattern,  respectively.  The  azimuthal  plane  is  defined  by 
0  =  PO*’  (the  coordinate  system  is  in  Fig.  1). 

Thetamin.  Thetamax.  Ntheta  REAL.  REAL.  INTEGER 

The  minimum  and  maximum  angles  of  0  (in  deg)  in  the  elevation 
plane,  and  the  number  of  points  which  are  calculated  in  the 
elevation  pattern,  respectively. 

Phiel _ REAL 

The  (j)  angle  of  the  elevation  plane  (in  deg).  In  other  words,  the 
elevation  plane  is  defined  by  (j)  =  Phiel. 


ThetaO.  PhiO _ REAL 

The  values  of  the  0  and  ([>  coordinates  (in  deg)  of  the  main  beam, 
respectively.  The  E-field  values  of  the  radiation  pattern  are 
normalized  by  the  amplitude  of  the  E-field  with  ThetaO  and  PhiO 
coordinates.  Furthermore,  ThetaO  and  PhiO  values  are  needed 
when  Nmode  =  4  (see  definition  of  the  Nmode),  and  they  are  used 
for  determining  the  phase  correction  of  the  current  on  each  patch. 
By  using  this  phase  correction  the  fields  from  each  patch  are  in 
phase  in  the  {ThetaO,  PhiO)  direction. 


Tvpepolar _ CHARACTER 

Typepolar  denotes  the  type  of  polarization  of  radiation  pattern.  L 
or  I  means  linear  polarization,  C  or  c  means  circular  polarization. 

Nmode _ INTEGER 

Nmode  is  the  mode  of  calculating  the  radiation  pattern.  Nmode=0,l 
and  5  means  that  the  radiation  pattern  is  calculated  rigorously,  i.e. 
the  results  from  the  SMiPA  program  are  used  (the  calculated 
current  distribution  is  used  for  chosen  frequency  point,  see 
explanation  of  Fmiri).  Nmode=0  means  calculation  of  the  gain 
pattern,  Nmode=l  means  that  the  radiation  pattern  is  normalized  to 
the  reference  direction,  and  Nmode=5  means  calculation  of  the 
element  pattern.  Nmode  2,  3  and  4  means  that  the  radiation  pattern 
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is  calculated  in  a  simple  way,  i.e.,  mutual  coupling  is  not  taken  into 
account  in  calculations  of  the  radiation  pattern.  When  Nmode  =  2  it 
is  supposed  that  patch  currents  have  equal  amplitudes  and  phases 
on  each  patch.  \^en  Nmode  =  3  amplitudes  and  phases  of  each 
patch  current  are  read  from  the  input  file  ‘SMiPApatin’.  When 
Nmode  =  4  all  patches  have  the  same  amplitude  of  the  current.  The 
phase  is  corrected  in  the  way  that  contributions  from  all  patches 
are  in  phase  in  the  direction  of  the  main  beam,  i.e.,  in  the  (thetaO, 
phiO)  direction.  If  Nmode=5,  then  the  number  of  the  patch  in  the  (j)- 
and  z-directions,  as  well  as  the  port  number  should  be  also  given  in 
the  same  line.  In  that  case  we  need  to  write  additional  three  integer 
numbers-coordinates). 


The  last  part  of  the  input  file  is  used  in  three  cases;  if  the  array  grid  is  neither  icosahedron  nor 
rectangular,  if  the  ports  do 'not  have  equal  excitations  and  we  want  to  run  SMiPA  program,  and  if 
we  want  to  calculate  the  radiation  pattern  (SMiPApat  program)  by  using  the  simple  method  and 
by  selecting  the  amplitudes  and  phases  of  currents  on  each  patch  {Nmode  =  3).  In  all  cases  the 
amplitudes  and  phases  of  the  excitation  currents,  as  well  as  patch  positions,  are  read  from  the 
end  of  the  input  file  SMiPA.in.  The  number  of  lines  in  the  file  is  the  same  as  the  number  of 
patches  times  number  of  ports.  Each  line  has  5  numbers:  first  two  numbers  are  the  patch  number- 
coordinates,  the  third  one  is  the  number  of  the  port  of  the  considered  patch,  and  the  fourth  and 
the  fifth  number  in  each  line  are  the  amplitude  and  the  phase  of  the  excitation  connected  to  the 
considered  port  of  the  patch. 

For  rectangular  grid  the  patch  number-coordinates  are  the  number  of  the  patch  in  0-  and  (]) 
direction  (the  patch  number  of  the  lowest  most  left  patch  is  (1,1)).  For  icosahedron  grid  the  patch 
number-coordinates  are  the  number  of  the  circle  and  the  number  of  the  patch  along  the  circle. 
For  non-rectangular  grid  first  two  numbers  represent  the  patch  position  (0  and  <|)  coordinate  in 
degrees).  We  can  check  if  these  last  lines  are  correctly  given  by  checking  the  description  of  the 
array  given  in  the  output  files  SMiPA.out  and  SMiPApat.out.  If  we  consider  some  other  case 
(e.g.  if  the  ports  have  equal  excitations)  there  is  no  need  for  these  input  data.  Therefore: 

If  needed,  next  Npat  x  Nport  lines  after  one  blank  line  are: 

Ipatl.  Ipat2.  Iport.  Amp,  Phase _ INTEGER.  REAL 

Ipatphi,  Ipatz,  Iport  are  integers  which  represents  the 
coordinates  of  each  port  (coordinate  Iport)  connected  to 
each  patch  (coordinates  Ipatl,  Ipat2).  For  that  port  real 
numbers  Amp  and  Phase  denote  the  amplitude  and  the 
phase  of  the  excitation  current. 


or 

Thetapat.  Phipat.  Iport.  Amp.  Phase  REAL.  INTEGER.  REAL 

Thetapat  and  Phipat  represent  the  0  and  <|)  coordinate  (in 
degrees)  of  the  center  of  each  patch,  Iport  is  the  number  of 
the  patch  port,  and  Amp  and  Phase  denote  the  amplitude 
and  the  phase  of  the  excitation  current. 
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5.1.3.2  Description  of  the  output  files 


There  are  six  output  files:  ‘SMiPA.out’,  ‘SMiPApat.ouf ,  ‘Sij.dat’,  ‘port.dat’  ‘pat-az.daf  and 
‘pat-el.daf ,  and  in  them  we  can  find  data  about  the  geometry  as  well  as  the  calculated  values  of 
the  input  impedance,  S-parameters  and  radiation  pattern  in  the  elevation  and  azimuthal  planes. 

SMiPA.out  The  file  contains  data  about  the  input  impedance  of  a  single  patch  (S  parameters  for 
array  case  are  given  in  Sij.dat).  When  constructing  the  array  we  need  first  to  determine  the 
parameters  of  a  single  patch,  e.g.  input  impedance,  return  loss,  and  impedance  bandwidth.  The 
structure  of  the  SMiPA  program  (and  output  file  ‘SMiPA.out’)  enables  these  calculations.  First 
part  of  the  file  contains  the  values  of  the  input  file.  Therefore  the  first  part  looks  similar  to  the 
input  file  ‘SMiPA.in’.  The  second  part  contains  the  input  impedance  and  return  loss  for  all 
calculated  frequencies.  Each  line  is  for  one  frequency  point,  and  it  has  four  values:  frequency  (in 
GHz),  normalized  input  impedance  (real  and  imaginary  part)  and  return  loss  (in  dB).  The 
calculated  input  impedance  is  normalized  by  the  characteristic  impedance  of  the  feeding 
transmission  line  ZO  (the  value  of  ZO  is  given  in  the  file  ‘SMiPA.in’;  for  aperture-coupled  patch 
antennas  it  is  calculated  from  the  geometry  of  feeding  structure). 

SMiPApatout  The  file  contains  data  about  the  array  geometry.  The  first  part  of  the  file 
contains  the  values  of  the  input  file.  Therefore  the  first  part  looks  similar  to  the  input  file 
‘SMiPA.in’.  The  second  part  describes  the  array  that  is  analyzed.  It  is  a  list  of  the  all  patch 
elements  in  the  array  with  their  coordinates  and  amplitude  and  phase  of  the  patch  current.  First 
two  columns  are  the  patch  number-coordinates,  the  third  column  is  the  number  of  the  basis 
function  on  the  considered  patch,  and  the  forth  and  fifth  columns  are  the  0-  and  <|)-coordinate  of 
the  patch.  The  sixth  and  seventh  columns  are  the  amplitude  and  phase  (in  deg)  of  the  considered 
basis  function,  respectively. 

Sij.dat  The  file  contains  the  calculated  S-parameters.  The  file  is  in  the  Touchstone  format,  i.e., 
using  this  file  it  is  possible  to  connect  SMiPA  and  Touchstone  programs.  The  first  line  is  a 
specification  line:  ‘#  GHz  S  DB  R’  and  at  the  end  of  the  line  there  is  a  value  of  the  normalization 
impedance.  GHz  means  that  the  frequency  values  are  in  GHz,  DB  means  that  S-parameters  are 
written  in  the  dB-angle  format.  Other  lines  in  the  file  are  the  S-parameters  for  different 
frequencies.  If  we  have  one  port  (one  patch)  then  each  line  contains  frequency,  |Sii|  in  dB,  and 
angle  of  Sn,  respectively.  In  the  case  of  two  ports  (two  patches)  each  line  contains  frequency, 
|Sn|  in  dB,  angle  of  Sn,  IS21I  in  dB,  angle  of  S21,  ISnl  in  dB,  angle  of  S12,  IS22I  in  dB,  and  angle  of 
S22,  respectively.  In  the  case  of  N  ports  each  frequency  is  described  with  N  lines,  each  containing 
|Sii|  in  dB,  angle  of  Sii,  |Si2|  in  dB,  angle  of  Si2,  etc.  Frequency  is  written  at  the  beginning  of  the 
first  of  A/^  lines. 

port.dat  The  file  contains  the  calculated  port  impedances  for  each  port.  For  each  frequency  point 
port  impedances  are  given  in  the  form  of  real  and  imaginary  part.  The  values  of  port  impedances 
are  normalized,  i.e.  port  impedances  are  divided  by  ZO. 

pat-az.dat  The  file  contains  the  field  values  in  the  azimuthal  plane  (the  azimuthal  plane  is 
defined  by  0=90*^).  The  input  values  of  variables  Phimin,  Phimax  and  Nphi  determine  the  ^ 
values  in  the  azimuthal  plane  for  which  the  radiation  pattern  is  calculated.  The  field  values  are 
normalized  with  the  field  value  with  ThetaO  and  PhiO  coordinates  that  are  defined  in  the  input 
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file  (the  exception  is  for  Nmode  =  0  and  Ninode  =  5).  The  first  and  second  columns  are  the  0  and 
(|)  coordinates  (in  deg)  for  which  the  radiation  pattern  is  calculated.  If  the  chosen  type  of 
polarization  is  linear  (parameter  Typepolar  in  the  input  file),  then  the  third  and  the  fourth 
columns  are  the  normalized  0-  and  (j)  -components  of  the  electric  field,  respectively.  If  the  chosen 
type  of  polarization  is  circular,  then  the  third  and  the  fourth  columns  are  the  normalized  right- 
hand  circular  polarization  (RHCP)  and  left-hand  circular  polarization  (LHCP)  components  of  the 
electric  field,  respectively. 

pat-el.dat  The  file  contains  the  field  values  in  the  elevation  plane  which  is  defined  by  ^-Phiel 
(the  value  of  Phiel  is  defined  in  the  input  file).  The  input  values  of  variables  Thetamin,  Thetamax 
and  Ntheta  determine  the  0  values  in  the  elevation  plane  for  which  the  radiation  pattern  is 
calculated.  The  field  values  are  normalized  with  the  field  value  with  coordinates  ThetaO  and 
PhiO  that  are  defined  in  the  input  file  (the  exception  is  for  Nmode  =  0  and  Nmode  =  5).  The  first 
and  second  columns  are  the  0  and  (|>  coordinates  (in  deg)  for  which  the  radiation  pattern  is 
calculated.  If  the  chosen  type  of  polarization  is  linear  (parameter  Typepolar  in  the  input  file), 
then  the  third  and  the  fourth  columns  are  the  normalized  0-  and  (j)  -components  of  the  electric 
field,  respectively.  If  the  chosen  type  of  polarization  is  circular,  then  the  third  and  the  fourth 
columns  are  the  normalized  right-hand  circular  polarization  (RHCP)  and  left-hand  circular 
polarization  (LHCP)  components  of  the  electric  field,  respectively. 
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5.2  Upgrade  of  the  Program  CyMPA 


5.2.1  Introduction 

This  chapter  gives  the  short  description  and  the  realized  improvements  for  of  the  program 
CyMPA.  The  abbreviation  CyMPA  means  the  program  for  analyzing  Cylindrical  Microstrip 
Patch  Arrays,  and  previous  version  of  the  program  was  designed  to  analyze  two  types  of  feeding 
structures:  coax  feeding  and  microstrip  line  feeding  [13].  The  upgraded  version  of  the  program 
can  also  analyze  aperture-coupled  cylindrical  patch  antennas.  Since  there  are  different 
requirements  on  calculation  of  radiation  pattern  and  input  impedance/mutual  coupling,  we  have 
actually  developed  two  programs  CyMPA  and  CyMPApat.  The  main  reason  for  splitting  the 
problem  into  two  parts  is  that  radiation  pattern  is  usually  calculated  only  for  one  frequency, 
while  input  impedance  and  mutual  coupling  is  usually  calculated  for  frequency  band  of  interest 
(and  thus  even  hundreds  of  frequency  points  can  be  calculated).  Therefore,  program  CyMPA 
calculates: 


•  current  distribution  at  each  patch  in  the  array  and  field  distribution  in  each  aperture 

•  input  impedance  at  each  input  port  in  the  array 

•  mutual  coupling  between  each  two  patches  in  the  array. 


and  program  CyMPApat  calculates 

•  radiation  pattern  of  the  array  when  all  mutual  couplings  are  taken  into  account 

•  radiation  pattern  of  the  array  without  taking  mutual  coupling  into  account  (fast  mode). 


In  order  to  determine  the  radiation  pattern,  program  CyMPApat  is  using  the  results  of  program 
CyMPA,  i.e.  current  distribution  on  patches  and  field  distribution  in  apertures  calculated  at 
selected  frequency  point  is  used. 

Program  is  written  in  FORTRAN  90,  and  therefore  can  be  compiled  and  run  at  each  machine 
where  FORTRAN  90  is  installed  or  which  support  the  executable  version  of  the  program.  The 
communication  with  the  program  is  made  via  input/output  ASCII  files.  For  case  of  simplicity, 
both  programs  CyMPA  and  CyMPApat  are  using  the  same  input  file  (CyMPA.in).  Furthermore, 
there  is  a  big  similarity  between  input/output  files  for  program  CyMPA  and  for  program  SMiPA. 
Therefore,  only  the  differences  between  cylindrical  and  spherical  case  will  be  explained  in  this 
chapter. 
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5.2.2  Differences  between  Cylindrical  and  Spherical  Program 

The  CyMPA  and  CyMPApat  programs  analyze  rectangular  patch  antennas  or  antenna  arrays 
placed  on  circular  cylinders,  as  illustrated  in  Fig.  5.3.  If  an  array  is  analyzed,  following 
assumptions  are  made: 

•  all  antenna  elements  are  identical 

•  the  patch  elements  are  of  rectangular  shape 

•  constant  spacing  between  neighboring  patches  in  circumferential  (^)  direction 

•  constant  spacing  between  neighboring  patches  in  axial  (z)  direction 

•  the  dielectric  structure  can  be  multilayered  with  losses 

•  the  patch  array  must  be  placed  at  the  top  of  one  such  layer 

•  in  the  case  of  aperture-coupled  patch  antennas  the  apertures  are  of  rectangular  shape  (with 
unequal  sides) 

•  there  must  be  a  constant  spacing  between  the  patch  center  and  the  aperture  center 


Figure  5.3.  The  geometry  of  the  cylindrical  patch  array 
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The  differences  between  SMiPA  and  CyMPA  and  their  input/output  files  are: 

•  In  the  cylindrical  case  the  patch  elements  are  of  rectangular  shape,  while  in  the  spherical  case 
patch  elements  can  be  both  rectangular  and  circular. 

•  In  the  cylindrical  case  only  the  rectangular  grid  of  patch  elements  can  be  considered.  In  the 
spherical  case  icosahedron  grid  can  be  also  considered,  as  well  as  arbitrary  grid. 

•  The  input  files  for  programs  CyMPA  and  SMiPA  are  called  cympa.in  and  smipa.in, 
respectively.  However,  their  structure  and  meaning  of  variables  is  very  similar. 

•  In  the  cylindrical  case  the  patch  width  /  feed  point  /  slot  width  /  grid  dimensions  are  given 
first  in  (|)-direction  and  then  in  z-direction.  In  the  spherical  case  dimensions  are  given  first  in 
0-direction  and  then  in  (|)-direction. 

•  Two  different  output  files  (but  with  almost  identical  meaning)  are  cympa.out  /  cympapat.out, 
and  smipa.out  /  smipapat.out.  All  other  output  files  have  the  same  meaning  and  structure. 


More  details  about  program  CyMPA  can  be  found  in  [13]. 
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6  BENCHMARKS 
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6.1  Spherical  16  element  Array  Following  the  Icosahedron  Grid 

-Radiation  by  Simple  Model 
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6.1.1  Introduction 

Icosahedron  grid  is  a  very  suitable  concept  for  spherical  arrays  that  steer  the  beam  throughout 
the  whole  hemisphere.  This  benchmark  presents  the  radiaition  properties  of  16  element 
microstrip  array  on  the  spherical  structure  following  icosahedric  grid,  as  depicted  in  figure  6. 1 . 
This  number  of  elements  corresponds  to  hemisphere  coverage  of  30  degrees.  The  step  of  15®  for 
each  next  circle  comes  from  dividing  by  4  the  side  of  triangle  of  basic  icosahedron  structure. 

The  calculations  are  performed  using  the  grounded  shell  of  radius  20  cm,  dimensions  of  patches 
are  2.5  cm  x  2.5  cm.  The  patch  excitation  for  each  patch  is  phase  corrected  in  order  to 
compensate  the  propagation  path  difference.  In  this  benchmark  input  file  and  calculations  for 
simple  and  fast  equiamplitude  model  are  presented.  This  model  gives  only  radiation  pattern  and 
does  not  include  mutual  coupling.  Radiation  pattern  calculation  procedure  when  using  this 
simple  approach  will  follow  the  diagram  in  Fig.  6.2.  Notice  that  only  SMiPApat  program  must 
be  applied.  The  next  benchmark  deals  with  the  same  configuration  calculated  rigorously,  so  that 
a  comparison  can  be  performed. 


Figure  6.1.  Spherical  array  with  16  elements  following  the  icosahedron  grid 


SMiPApat .  out 


OUTPUT  DATA 


RADIATION  PATTERN 
IN  ELEVATION  PLANE 


Figure  6.2.  Procedure  for  calculating  radiation  pattern  using  simple  approach 
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6.1.2  Summary  technical  description 


Antenna  /  Antenna  Array  description 


Type:  16  element  spherical  array  following  icosahedron  grid 


Frequency:  3.64  GHz 


The  mechanical  and  dielectric  data  of  the  structure: 


no,  of 
layer 

layer 

description 

relative 
permi t ti vi ty 

loss 

tcuigent 

inner 

radius 

outer 

radius 

thickness 

- 

ground 

sphere 

- 

- 

20  cm 

20  cm 

- 

1 

PTFE  sphere 

2.52 

0.0019 

20  cm 

20.1576  cm 

0.1576  cm 

Rectangular  patch  element  data: 


Element  is  printed  on 

azimuthal 

dimension 

elevation 

dimension 

outer  shell  of 
layer  no .  1 

2.5  cm 

2.5  cm 

Coaxial  excitation  position,  relative  to  patch  center: 


elevational  position 

azimuthal  position 

0 . 3  cm 

0  cm 

Array  configuration:  16  element  icosahedron  grid 


number  of  rings 

icosahedron  side 
subdivision  factor 

3 

4 
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6.1.3  Numerical  results 

The  radiation  pattern  in  E  and  H  plane  were  calculated  using  assumption  of  equiamplitude 
identical  patch  currents  that  does  not  take  mutual  coupling  into  account.  The  calculation  is 
performed  at  frequency  of  3.64  GHz.  Figure  6.3  shows  the  input  file  for  this  case.  Light  printed 
are  the  values  and  rows  that  are  irrelevant  for  this  case. 

The  SMiPApat  routine  was  run  to  calculate  the  radiation  pattern.  The  calculated  radiation  pattern 
in  E  and  H  plane  is  given  in  Figure  6.4.  and  Figure  6.5. 


3.64 

20.0 

1 

3.64 

1  R 

?  Frequency  (GHz) 

!  Ground  shell  and  patch  radius  (cm) 

!  Number  of  dielectric  layers 

0.1576 

1 

2.52 

0.0019 

!  Relative  permittivity  and  loss  tanges 
!  Die.  layer  at  which  patches  are  placed 

3 

4 

I 

!  Array  grid  parameters  and  grid  type 

2.5 

C 

1 

2.5 

R 

!  Patch  width  in  theta  and  phi  dir (cm) ,rectang. sh 
!  Feeding  type:  'cV=cocUc,  ’m'=microstrip,  *a*=aportu 
!  Number  of  ports  per  patch 

0.3 

Y 

0.0 

!  Port  position 

?  Equal  port  excitations?  (Y  -  yes^  N  -  no) 

0.0 

50.0 

0.0 

!  Distance  between  patch  centers 
{  Characteristic  in^edance  of  the  feed  line 

-180.0 

180.0 

361 

!  Azimuthal  pattern 

-89.0 

0.0 

269.0 

359 

!  Elevation  pattern 

!  Phi  angle  of  the  evelation  pattern 

90.0 

L 

4 

0.0 

!  Theta  and  phi  angle  of  the  main  beam 
!  Polarization  (L  -  linear,  C  -  circular) 

!  Phase  corrected  excitation 

Figure  6.3.  Input  file  smipa.in 


A  few  lines  from  the  radiation  pattern  output  file  ‘pat-az.dat’  around  (|)  =  0  are: 


90.0000 

-5.0000 

90.0000 

-4.0000 

90.0000 

-3.0000 

90.0000 

-2.0000 

90.0000 

-1.0000 

90.0000 

0.0000 

90.0000 

1.0000 

90.0000 

2.0000 

90.0000 

3.0000 

90.0000 

4.0000 

90.0000 

5.0000 

-0.7020 

-200.0000 

-0.4477 

-200.0000 

-0.2512 

-200.0000 

-0.1114 

-200.0000 

-0.0278 

-200.0000 

0.0000 

-200.0000 

-0.0278 

-200.0000 

-0.1114 

-200.0000 

-0.2512 

-200.0000 

-0.4478 

-200.0000 

-0.7021 

-200.0000 

relative  power  [dB]  relative  power  [dB] 
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azimuthal  angle  [deg] 

Figure  6.4.  Calculated  radiation  pattern  in  azimuthal  plane 


elevation  angle  [deg] 


Figure  6.5.  Calculated  radiation  pattern  in  elevation  plane 
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6.2  Spherical  16  element  Array  Following  the  Icosahedron  Grid 

-Radiation  by  Rigorous  Model 
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6.2.1  Introduction 

In  the  previous  benchmark  we  considered  the  radiation  properties  of  16  element  microstrip  array 
on  the  spherical  structure  following  icosahedric  grid  (Fig.  6.6).  The  radiation  pattern  was 
calculated  without  taking  mutual  coupling  into  account.  In  this  benchmark  the  radiation  pattern 
calculated  in  a  rigorous  way  will  be  considered.  Therefore,  after  running  SMiPA  program  and 
obtaining  currents  on  each  patch,  SMiPApat  program  must  be  run  in  order  to  obtain  radiation 
pattern.  The  comparison  of  radiation  pattern  with  the  simple  calculation  model  that  is  used  in 
previous  benchmark  is  performed. 

The  calculations  are  performed  using  the  grounded  shell  of  radius  20  cm,  dimensions  of  patches 
are  2.5  cm  x  2.5  cm.  The  patch  excitation  for  each  patch  is  phase  corrected  in  order  to 
compensate  the  propagation  path  difference.  16  elements  correspond  to  hemisphere  coverage  of 
30  degrees.  The  step  of  15°  for  each  next  circle  comes  from  dividing  by  4  the  side  of  triangle  of 
basic  icosahedron  structure. 

Radiation  pattern  calculation  procedure,  when  using  rigorous  approach,  will  follow  the 
diagram  in  Fig.  6.7.  SMiPA  program  must  be  applied  to  determine  rigorous  current  distributions 
on  the  patch  before  calculating  radiation  pattern.  In  order  to  obtain  radiation  pattern  SMiPApat 
program  must  be  run. 


Figure  6.6.  Spherical  array  with  16  elements  following  the  icosahedron  grid 
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SMiPA.out 
INPUT  DATA  I 

sij .dat 
S  PARAMETERS 


port . dat 


Figure  6.7.  Procedure  for  calculating  radiation  pattern  using  rigorous  approach 
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6.2.2  Summary  technical  description 


Antenna  /  Antenna  Array  description 


Type:  16  element  spherical  array  following  icosahedron  grid 


Frequency:  3.64  GHz 


The  mechanical  and  dielectric  data  of  the  structure: 


no.  of 
layer 

layer 

description 

relative 
permi t ti vi ty 

loss 

tangent 

inner 

radius 

outer 

radius 

thickness 

- 

ground 

sphere 

- 

- 

20  cm 

20  cm 

- 

1 

PTFE  sphere 

2.52 

0.0019 

20  cm 

20.1576  cm 

0.1576  cm 

Rectangular  patch  element  data: 


Element  is  printed  on 

azimuthal 

dimension 

elevation 

dimension 

outer  shell  of 
layer  no .  1 

2.5  cm 

2.5  cm 

Coaxial  excitation  position,  relative  to  patch  center: 


elevational  position 

azimuthal  position 

0.3  cm 

0  cm 

Array  configuration:  16  element  icosahedron  grid 


number  of  rings 

icosahedron  side 
subdivision  factor 

3 

4 
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6.2.3  Numerical  results 

The  radiation  pattern  in  E  and  H  plane  were  calculated  using  assumption  of  equiamplitude 
identical  patch  currents  that  does  not  take  mutual  coupling  into  account.  The  calculation  is 
performed  at  frequency  of  3.64  GHz.  Figure  6.8  shows  the  input  file  for  this  case.  Light  printed 
are  the  values  and  rows  that  are  irrelevant  for  this  case. 

The  SMiPApat  routine  was  run  to  calculate  the  radiation  pattern.  The  calculated  radiation  pattern 
in  E  and  H  plane  is  given  in  Figure  6.9. 


3.64 

3.64 

1 

R  ! 

Frequency  (GHz) 

20.0 

f 

Ground  shell  and  patch  radius  (cm) 

1 

! 

Number  of  dielectric  layers 

0.1576 

2.52 

0.0019  ! 

Relative  permittivity  and  loss  tcuiges 

1 

t 

Die.  layer  at  which  patches  are  placed 

3 

4 

I 

! 

Array  grid  parameters  and  grid  type 

2.5 

2.5 

R 

1 

Patch  width  in  theta  and  phi  dir  (cm)  ,recteuig.sh 

C 

! 

Feeding  type : *  c  *  =coax , ' m ' =micros trip , *  a ' =apertu 

1 

! 

Number  of  ports  per  patch 

0.3 

0.0 

1 

Port  position 

N 

1 

Equal  port  excitations?  (Y  -  yes,  N  -  no) 

0.0 

0.0 

1 

Distance  between  patch  centers 

50.0 

( 

Characteristic  impedance  of  the  feed  line 

-180.0 

180.0 

361 

1 

Azimuthal  pattern 

-89.0 

269.0 

359 

1 

Elevation  pattern 

0.0 

f 

Phi  angle  of  the  evelation  pattern 

90.0 

0.0 

f 

Theta  and  phi  euigle  of  the  main  becun 

L 

f 

Polarization  (L  -  linear,  C  -  circular) 

1 

! 

Radiation  pattern  calculated  rigorously 

1  1 

1 

1.0 

0.0 

2  1 

1 

1.0 

30.0 

2  2 

1 

1.0 

30.0 

2  3 

1 

1.0 

30.0 

2  4 

1 

1.0 

30.0 

2  5 

1 

1.0 

30.0 

3  1 

1 

1.0 

118.0 

3  2 

1 

1.0 

118.0 

3  3 

1 

1.0 

118.0 

3  4 

1 

1.0 

118.0 

3  5 

1 

1.0 

118.0 

3  6 

1 

1.0 

118.0 

3  7 

1 

1.0 

118.0 

3  8 

1 

1.0 

118.0 

3  9 

1 

1.0 

118.0 

3  10 

1 

1.0 

118.0 

Figure  6.8.  Input  file  smipa.in 


A  few  lines  from  the  radiation  pattern  output  file  'pat-az.dat’  around  ^  =  0  are: 


90.0000 

90.0000 

90.0000 

90.0000 

90.0000 


-5.0000 

-4.0000 

-3.0000 

-2.0000 

-1.0000 


-0.7268 

-0.4627 

-0.2588 

-0.1142 

-0.0281 


-79.3460 

-81.1226 

-83.3903 

-86.4056 

-90.4843 
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90.0000 

0.0000 

0.0000 

-93.8173 

90.0000 

1.0000 

-0.0299 

-90.8304 

90.0000 

2.0000 

-0.1179 

-86.8503 

90.0000 

3.0000 

-0.2645 

-83.9390 

90.0000 

4.0000 

-0.4705 

-81.7928 

90.0000 

5.0000 

-0.7369 

-80.1506 

azimuthal  angle  [deg] 


(a) 


elevation  angle  [deg] 

(b) 

Figure  6.9.  Comparison  of  calculated  radiation  pattern  in  elevation  plane  using  simple  and 
rigorous  models,  (a)  azimuthal  plane,  (b)  elevation  plane. 
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6.3  Rectangular  patch  antenna  on  spherical  structure 
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6.3.1  Introduction 

This  example  deals  with  input  impedance  of  rectangular  patch  printed  on  spherical  structure.  The 
geometry  of  the  problem  is  given  in  Fig.  6.10.  The  square  patch  of  dimension  5.1  x  5.1  cm  is 
placed  above  the  ground  shell  using  syrofoam  spacers.  Therefore,  substrate  relative  dielectric 
constant  is  8^  =  1 .0  and  the  air  gap  due  to  syrofoam  spacers  is  0.52  cm  thick.  For  the  validation 
of  calculation  results  a  laboratory  model  was  built  from  a  copper  half-sphere  of  radius  a=  1 8.7 
cm  at  which  patches  can  be  mounted  at  arbitrary  positions. 

The  detailed  data  sheet  about  the  calculated/measured  antenna  is  given  in  6.3.2.  In  order  to 
calculate  this  example  a  run  of  SMiPA  is  needed,  with  an  input  file  that  is  shown  in  section 
6.3.3.  The  comparison  of  calculation  and  measurement  can  be  found  in  section  6.3.3  as  well.  The 
calculated  and  measured  impedances  are  in  a  good  agreement. 


Figure  6. 10.  Rectangular  patch  antenna  on  spherical  structure 
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6.3.2  Summary  technical  description 


Antenna  /  Antenna  Array  description 


Type:  Rectangular  patch  element  printed  on  spherical  grounded  substrate 


Frequency:  2  GHz  -  3  GHz,  21  equidistant  points 


The  mechanical  and  dielectric  data  of  the  structure: 


no.  of 
layer 

layer 

description 

relative 

peinoittivity 

loss 

tcuigent 

inner 

radius 

outer 

radius 

thickness 

- 

ground 

sphere 

- 

18.7  cm 

18.7  cm 

- 

1 

air  gap 

o 

o 

0.0000 

18.7  cm 

19.22  cm 

0.52  cm 

Rectangular  patch  element  data: 


Element  is  printed  on 

azimuthal 

dimension 

elevation 

dimension 

outer  shell  of 
layer  no .  1 

5.1  cm 

5 . 1  cm 

Coaxial  excitation  position,  relative  to  patch  center: 


elevational  position 

azimuthal  position 

1.95  cm 

0  cm 

Array  configuration:  single  element 
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6.3.3  Numerical  and  experimental  results 

The  calculation  is  performed  in  frequency  range  of  2.0  to  3.0  GHz  at  21  equidistant  frequency 
points.  Figure  6.1 1  shows  the  input  file  for  the  program  SMiPA.  Light  printed  are  the  values  and 
rows  that  are  irrelevant  for  this  case.  Here,  distance  between  patches  is  irrelevant  since  a  single 
patch  configuration  is  considered.  Additionally,  all  data  for  radiation  pattern  is  irrelevant,  since 
in  this  benchmark  only  the  input  impedance  is  of  interest. 


2.0 

18.7 

1 

3.0 

21  R 

!  Frequency  (GHz) 

!  GroTjnd  shell  and  patch  radius  (cm) 

!  Number  of  dielectric  layers 

0.52 

1 

1.0 

0.0000 

!  Relative  permittivity  and  loss  tanges 
!  Die.  layer  at  which  patches  are  placed 

1 

1 

R 

!  Array  grid  parameters  and  grid  type 

5.1 

C 

1 

5.1 

R 

!  Patch  width  in  theta  and  phi  dir (cm) ,rectang.sh 
!  Feeding  type : *  c  *  =coax , ' m ’ =micro strip , ’ a ’ =apertu 
!  Number  of  ports  per  patch 

1.95 

Y 

o 

o 

!  Port  position 

!  Equal  port  excitations?  (Y  -  yes,  N  -  no) 

o  o 

o  o 
in 

o 

o 

!  Distance  between  patch  centers 

!  Characteristic  impedcince  of  the  feed  line 

-180.0 

180.0 

361 

\  Azimuthal  pattern 

-89.0 

0.0 

269.0 

359 

?  Elevation  pattern 

!  Phi  angle  of  the  evelation  pattern 

90.0 

I. 

1 

o 

o 

!  Theta  and  phi  angle  of  the  main  beam 
?  Polarization  (L  -  linear,  C  -  circular) 

!  Phase  corrected  excitation 

Figure  6.11.  Input  file  sm  ipa.  in 


A  few  lines  from  the  file  ‘SMiPA.out’  are: 


2.35000 

0.65205 

3.13193 

-1.01286 

2.40000 

1.01869 

3.43469 

-1.28849 

2.45000 

1.67610 

3.72735 

-1.66488 

2.50000 

2.80854 

3.72161 

-2.19101 

2.55000 

4.14682 

2.68146 

-2.94561 

2.60000 

4.23653 

0.69463 

-4.05946 

2.65000 

3.05476 

-0.43903 

-5.76081 

2.70000 

1.98459 

-0.58367 

-8.48768 

2.75000 

1.32226 

-0.39931 

-13.24033 

The  comparison  of  the  calculated  input  impedance  is  given  in  Fig.  6.12  for  the  real  part  and  in 
Fig.  13  for  the  imaginary  part.  The  measured  results  are  also  given. 


normalized  impedance  normalized  impedance 
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Figure  6.13.  Comparison  of  calculated  and  measured  imaginary  part  of  input  impedance 
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6.4  Circular  patch  antenna  on  spherical  structure 
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6.4.1  Introduction 

This  example  deals  with  input  impedance  of  circular  patch  printed  on  the  spherical  structure.  The 
geometry  of  the  problem  is  given  in  figure  6.14.  The  analyzed  circular  patch  has  radius  of  1.45 
cm,  and  it  is  printed  on  the  substrate  with  dielectric  constant  e  ,  =  2.52  and  thickness  h  -  1.576 
mm. 

The  detailed  data  sheet  about  the  calculated/measured  antenna  is  given  in  6.4.2.  In  order  to 
calculate  this  example  a  run  of  SMiPA  is  needed,  with  an  input  file  that  is  shown  in  section 
6.4.3.  The  influence  of  the  radius  of  the  spherical  structure  on  the  input  impedance,  cases  with 
several  different  radii  (2.5  cm,  4  cm  and  10  cm)  are  calculated  and  the  results  compared  with 
measured  results  of  their  planar  counterpart. 


Figure  6. 14.  Circular  patch  antennas  on  spherical  structure 
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6.4.2  Summary  technical  description 


Antenna  /  Antenna  Array  description 

Type:  Circular  patch  element  printed  on  spherical  grounded  substrate 
Frequency:  3.2  GHz -4.0  GHz,  81  equidistant  points 
The  mechanical  and  dielectric  data  of  the  structure: 


no.  of 
layer 

layer 

description 

relative 

permittivity 

loss 

tangent 

inner 

radius 

outer 

radius 

thickness 

ground 

sphere 

- 

- 

planar 

planar 

- 

1 

PTFE  sphere 

2.52 

0.0019 

planar 

planar 

0.1576  cm 

no.  of 
layer 

layer 

description 

relative 

permittivity 

loss 

tangent 

inner 

radius 

outer 

radius 

thickness 

- 

ground 

sphere 

- 

- 

10  cm 

10  cm 

- 

1 

PTFE  sphere 

2.52 

0.0019 

10  cm 

10.1576  cm 

0.1576  cm 

no.  of 
layer 

layer 

description 

relative 
permi t t i vi ty 

loss 

tangent 

outer 

radius 

thickness 

- 

ground 

sphere 

- 

- 

4  cm 

4  cm 

- 

1 

PTFE  sphere 

2.52 

0.0019 

4  cm 

4.1576  cm 

0.1576  cm 

no.  of 
layer 

layer 

description 

relative 
permi  tti  vi  ty 

loss 

tangent 

outer 

radius 

thickness 

j 

ground 

sphere 

- 

- 

2.5  cm 

2.5  cm 

1 

PTFE  sphere 

2.52 

0.0019 

2.5  cm 

2.6576  cm 

0.1576  cm 

Circular  patch  element  data: 


Element  is  printed  on 

diameter 

outer  shell  of 
layer  no .  1 

2.9  cm 

Coaxial  excitation  position,  relative  to  patch  center:  0.3  cm 


Array  configuration:  single  element 
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6.4.3  Numerical  and  experimental  results 

The  calculation  is  performed  in  frequency  range  of  3.2  to  4.0  GHz  at  81  equidistant  frequency 
points.  Figure  6.15  shows  the  input  file  for  the  program  SMiPA  in  case  of  10  cm  diameter.  Light 
printed  are  the  values  and  rows  that  are  irrelevant  for  this  case.  Here,  distance  between  patches  is 
irrelevant  since  a  single  patch  configuration  is  considered.  Additionally,  all  data  for  radiation 
pattern  is  irrelevant,  since  in  this  benchmark  only  the  input  impedance  is  of  interest. 


3.2 

4.0 

81 

R  ?  Frequency  (GHz) 

10.0 

!  Groxand  shell  and  patch  radius  (cm) 

1 

!  Number  of  dielectric  layers 

0.1576 

2.52 

0.0019 

!  Relative  permittivity  and  loss  tanges 

1 

!  Die.  layer  at  which  patches  are  placed 

1 

1 

R 

!  Array  grid  parameters  and  grid  typo 

2.9 

4.0 

C 

!  Patch  diameter (cm) , circular  shape 

C 

!  Feeding  type : ' c  * -coax , ' m ' =micros trip , ’ a ' =apertu 

1 

!  Number  of  ports  per  patch 

0.3 

o 

o 

!  Port  position 

Y 

!  Equal  port  excitations?  (Y  -  yes,  N  -  no) 

0.0 

0.0 

t  Distance  between  patch  centers 

50.0 

!  Characteristic  impedance  of  the  feed  line 

*~180.0 

180.0 

361 

!  Azimuthal  pattern 

-89.0 

269,0 

359 

!  Elevation  pattern 

0.0 

!  Phi  angle  of  the  evelation  pattern 

o 

o 

<3\ 

0.0 

!  Theta  and  phi  angle  of  the  main  beam 

L 

1  Polarization  (L  -  linear,  C  -  circular) 

1 

1  Phase  corrected  excitation 

Figure  6. 15.  Input  file  smipa.m  for  10  cm  ground  shell  diameter 


A  few  lines  from  the  file  ‘smipa.out’  are: 


3.62000 

0.37823 

0.88159 

-3.61726 

3.63000 

0.45047 

0.89057 

-4.22490 

3.64000 

0.53641 

0.88651 

-4.97472 

3.65000 

0.63469 

0.86152 

-5.90979 

3.66000 

0.73924 

0.80634 

-7.09043 

3.67000 

0.83714 

0.71356 

-8.60387 

3.68000 

0.90922 

0.58359 

-10.57906 

3.69000 

0.93675 

0.42977 

-13.19239 

3.70000 

0.91200 

0.27617 

-16.47603 

3.71000 

0.84339 

0.14638 

-18.71635 

3.72000 

0.74956 

0.05297 

-16.69852 

3.73000 

0.64904 

-0.00401 

-13.43902 

3.74000 

0.55414 

-0.03178 

-10.82575 

3.75000 

0.47065 

-0.03914 

-8.85473 

3.76000 

0.39997 

-0.03357 

-7.34775 

3.77000 

0.34128 

-0.02046 

-6.17312 

The  comparison  of  the  calculated  input  impedances  is  given  in  Fig.  6.16  for  the  real  part  and  in 
Fig.  17  for  the  imaginary  part.  The  measured  results  for  planar  case  are  also  given. 
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r  =  2.5cm 
r  =  4  cm 
r =  10cm 
planar  (meas.) 


frequency  (GHz) 

Figure  6.16.  Comparison  of  calculated  real  part  of  input  impedance  for  spheres  with  radii  of 
2.5,  4  and  10  cm  with  measurements  for  the  planar  case 


E 

T  O.B 


i  0.2 


—  r  =  2.5cm 

—  -  r  =  4  cm 

—  r=10cm 
-  planar  (meas.) 


\ 


'3.2  3.3  3.4  3.5  3.B  3.7  3.8  3.9  4 

frequency  (GHz) 


Figure  6. 1 7.  Comparison  of  calculated  imaginary  part  of  input  impedance  for  spheres  with 
radii  of  2.5,  4  and  10  cm  with  measurements  for  the  planar  case 


6.5  Mutual  Coupling  of  Two  Rectangular  Patches 

Printed  on  Spherical  Structure 
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6.5.1  Introduction 

This  example  deals  with  mutual  coupling  of  two  rectangular  patches  printed  on  spherical 
structure.  The  geometry  of  the  problem  is  given  in  Fig.  6.18.  The  square  patches  of  dimension 

5.1  X  5.1  cm  are  placed  above  the  ground  shell  using  syrofoam  spacers.  Therefore  substrate 
relative  dielectric  constant  is  =  1.0  and  the  air  gap  due  to  syrofoam  spacers  is  0.52  cm  thick. 
For  the  validation  of  calculation  results  a  laboratory  model  was  built  from  a  copper  half-sphere 
of  radius  a  =  1 8.7  cm  at  which  patches  can  be  mounted  at  arbitrary  positions. 

The  detailed  data  sheet  about  the  calculated/measured  antenna  is  given  in  6.5.2.  In  order  to 
calculate  this  example  a  run  of  SMiPA  is  needed,  with  an  input  file  that  is  shown  in  section 
6.5.3.  The  comparison  of  calculation  and  measurement  can  be  found  in  section  6.3.3  as  well.  The 
calculated  and  measured  mutual  coupling  is  in  a  good  agreement. 


Figure  6. 18.  Two  element  array  of  rectangular  patches  on  spherical  structure 
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6.5.2  Summary  technical  description 


Antenna  /  Antenna  Array  description 


Type:  Rectangular  patch  element  printed  on  spherical  grounded  substrate 


Frequency:  2  GHz  -  3  GHz,  21  equidistant  points 


The  mechanical  and  dielectric  data  of  the  structure: 


no.  of 
layer 

layer 

description 

relative 
permi  t  ti  vi  ty 

loss 

tangent 

inner 

radius 

outer 

radius 

thickness 

- 

ground 

sphere 

- 

- 

18.7  cm 

18.7  cm 

-- 

1 

air  gap 

1.00 

0.0000 

18.7  cm 

19.22  cm 

0.52  cm 

Rectangular  patch  element  data: 


Element  is  printed  on 

azimuthal 

dimension 

elevation 

dimension 

outer  shell  of 
layer  no .  1 

5.1  cm 

5.1  cm 

Coaxial  excitation  position,  relative  to  patch  center: 


elevational  position 

azimuthal  position 

1 . 95  cm 

0  cm 

Array  configuration: 


patch  center  distances 

no. of  elements 

separation  direction 

10.25  cm 

2 

6  (azimuthal) 
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6.5.3  Numerical  results 

The  calculation  is  performed  in  frequency  range  of  2.0  to  3.0  GHz  at  21  equidistant  frequency 
points.  Figure  6. 19  shows  the  input  file  “smipa.in”  for  the  program  SMiPA.  Light  printed  are  the 
values  and  rows  that  are  irrelevant  for  this  case.  All  data  for  radiation  pattern  is  irrelevant,  since 
in  this  benchmark  only  the  input  impedances  and  mutual  coupling  are  of  interest. 


2.0 

18.7 

1 

3.0 

21  R 

!  Frequency  (GHz) 

!  Ground  shell  and  patch  radius  (cm) 

!  Number  of  dielectric  layers 

0.52 

1 

1.0 

0.0000 

!  Relative  permittivity  and  loss  tanges 
!  Die.  layer  at  which  patches  are  placed 

2 

1 

R 

!  Array  grid  parameters  cuid  grid  type 

5.1 

C 

1 

5.1 

R 

!  Patch  width  in  theta  and  phi  dir (cm) ,rectang. sh 
!  Feeding  type: 'c'=coax,  *m*=:micr  os  trip, 'a'=apertu 
!  Number  of  ports  per  patch 

1.95 

Y 

0.0 

!  Port  position 

!  Equal  port  excitations?  (Y  -  yes,  N  -  no) 

10.25 

50.0 

0.0 

!  Distance  between  patch  centers 
!  Characteristic  impedance  of  the  feed  line 

-180.0 

180.0 

361 

>  Azimuthal  pattern 

-89.0 

0.0 

269.0 

359 

!  Elevatxon  pattern 

!  Phi  angle  of  the  evelation  pattern 

90.0 

L 

1 

0.0 

!  Theta  and  phi  angle  of  the  main  beam 

1  Polarization  (L  -  linear,  C  -  circular) 

!  Phase  corrected  excitation 

Figure  6. 19.  Input  file  cimpa.in 


A  few  lines  from  the  file  ‘Slj.dat’  are: 


2.2500 

-0.65 

40.60 

-39.53 

146.36 

-39.53 

146.36 

-0.65 

40.59 

2.3000 

-0.80 

37.61 

-38.48 

141.07 

-38.48 

141.07 

-0.81 

37.60 

2.3500 

-1.01 

34.15 

-37.33 

135.41 

-37.33 

135.41 

-1.01 

34.15 

2.4000 

-1.28 

30.13 

-36.05 

129.18 

-36.05 

129.18 

-1.29 

30.13 

2.4500 

-1.66 

25.39 

-34 . 63 

122.07 

-34.63 

122.07 

-1.67 

25.39 

2.5000 

-2.18 

19.72 

-33.08 

113.68 

-33.08 

113.68 

-2.19 

19.74 

2.5500 

-2.94 

12.86 

-31.42 

103.54 

-31.42 

103.54 

-2.95 

12.90 

2.6000 

-4.05 

4.46 

-29.71 

91.13 

-29.71 

91.13 

-4.06 

4.52 

2.6500 

-5.75 

-6.05 

-28.06 

76.04 

-28.06 

76.04 

-5.77 

-5.94 

2.7000 

-8.49 

-19.89 

-26.61 

58.17 

-26.61 

58.17 

-8.50 

-19.70 

The  comparison  of  the  calculated  mutual  coupling  with  the  measurements  is  given  in  Fig.  6.20. 


Analysis  of  Curved  Microstrip  Antennas 


85 


Figure  6.20.  Calculated  and  measured  mutual  coupling  coefficient  S21. 
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6.6  Aperture-Coupled  Cylindrical  Patch  Antenna 
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6.6.1  Introduction 

The  developed  program  for  analyzing  cylindrical  aperture-coupled  patches  has  been  verified 
against  measurements  made  at  University  of  Karlsruhe  [16].  The  aperture-coupled  array  was 
realized  in  C-band  ,  and  it  consists  of  supporting  metallic  cylinder  with  radius  2.6  cm,  covered 
with  foam  spacer  of  thickness  1.45  mm  and  with  RT/duroid  substrate.  On  RT/duroid  substrate 
both  microctrip  transmission  line  (facing  supporting  metallic  cylinder  covered  with  foam)  and 
ground  plane  with  apertures  (facing  outwards)  were  printed.  The  outer  region  consists  of  foam 
spacer  (thickness  0.4  cm)  and  Kapton  foil  on  which  patches  were  printed.  Notice  that  the  radius 
of  the  ground  tube  with  apertures  is  calculated  as  a  sum  of  the  radius  of  the  supporting  metallic 
cylinder  and  of  the  thicknesses  of  the  foam  and  RT/duroid  substrates.  The  patch  size  is  2.5  cm 
(azimuthal  dim.)  x  1.99  cm  (axial  dim.),  the  slot  size  is  1.38  x  0.09  cm,  the  length  of  the  open 
stub  is  0.35  cm,  and  the  width  of  the  microstrip  line  is  0.154  cm.  The  2x2  array  is  considered 
and  as  a  verification  Sn  parameter,  S21  parameter  and  radiation  pattern  are  calculated  and 
compared  with  measurements.  The  patches  are  separated  in  (Ji-direction  with  distance  between 
neighboring  patch  centers  of  3.54  cm  and  in  0-direction  with  distance  between  neighboring  patch 
centers  of  3.46  cm.  The  design  under  analysis  is  shown  in  Fig.  6.21. 


Figure  6.21  -  The  circular-cylindrical  2x2  aperture-coupled  patch  array 
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6.6.2  Summary  technical  description 

Measurements:  measurements  made  at  University  of  Karlsruhe  [10] 

Antenna  /  Antenna  Array  description 
Type:  2x2  element  array 

Frequency:  C-band,  4.5  GHz  -  6.5  GHz,  calculation:  1 1  equally  spaced  points 
The  mechanical  and  dielectric  data  of  the  structure: 


no.  of 
layer 

layer 

description 

relative 

permittivity 

loss 

tangent 

inner 

radius 

outer 

radius 

thickness 

inner 

2 

foam 

spacers 

1.1 

0.001 

2.6  cm 

4.05  cm 

1.45  cm 

inner 

1 

RT/duroid 

siibstrate 

2.2 

0.0009 

4 . 05  cm 

4.1008  cm 

0.0508  cm 

- 

ground 

tube 

- 

4.1008  cm 

4 . 1008  cm 

- 

1 

foam 

spacers 

1.1 

0.001 

4.1008  cm 

4 . 5008  cm 

0.4  cm 

2 

Kapton 

foil 

3.7 

0.001 

4.5008  cm 

4.5083  cm 

0.0075  cm 

Rectangular  patch  element  data: 


Element  is  printed  on 

azimuthal 

dimension 

axial 

dimension 

outer  shell  of 
layer  no .  2 

2.5  cm 

1.99  cm 

Aperture  coupled  excitation  data: 


aperture 

(slot) 

length 

aperture 

(slot) 

width 

microstrip 
line  width 

open 

stub 

length 

substrate 

thickness 

Er 

sxibstrate 

loss 

tangent 

1.38  cm 

0.09  cm 

0.154  cm 

0.35  cm 

0.0508  cm 

2.2 

0.0009 

Array  configuration: 


patch  center  distance 

no. of  elements 

separation  direction 

3.54  cm 

2 

^  (azimuthal) 

3.46  cm 

2 

z  (axial) 
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6.6.3  Numerical  results 

As  a  verification  case  we  wanted  to  calculate  Sn  parameter,  S21  parameter  and  radiation  pattern 
of  a  2  X  2  patch  array.  The  input  file  “cimpa.in”  for  this  case  is  shown  in  figure  6.22. 


|4.5  6.5 

11  R 

!  Frequency  (GHz) 

4.1008 

!  Ground  tube  radius  (cm) 

2 

!  Number  of  dielectric  layer 

0.4 

1.100 

0.001 

!  h,  epsilon_r,  tan(delta) 

0.0075 

3.700 

0.001 

!  h,  epsilon_r,  tan(delta) 

2 

!  Die.  layer  at  which  patches  are  placed 

2 

2 

R 

!  Number  of  patches  in  phi/z  direction 

2.5 

1.99 

!  Patch  width  in  phi  and  z  direction  (cm) 

A 

!  Feeding  type:  ‘c*=coax,  *m’=microstrip,  ‘a’=aperture 

0.0 

0.0 

!  Aperture  position  in  phi/z  direction  (cm) 

Y 

!  Equal  port  excitations?  (Y  -  yes,  N  -  no) 

3.54081 

3.46154 

!  Distance  between  patch  centers 

1.38 

0.09 

!  Slot  dimensions 

2 

!  Number  of  dielectric  layers  inside 

0.0508 

2.2 

0.0009 

!  h,  epsilon_r,  tan(delta) 

4.05 

1.1 

0.001 

!  h,  epsilon_r,  tan(delta) 

1 

!  Die.  layer  at  which  feed  line  is  placed 

0.154 

0.35 

!  Feed  line  width  and  stub  length 

-180.0 

180.0 

361 

!  Azimuthal  pattern 

1.0 

179.0 

179 

!  Elevation  pattern 

0.0 

!  Phi  angle  of  the  evelation  pattern 

90.0 

0.0 

!  Theta  and  phi  angle  of  the  main  beam 

L 

!  Polarization  (L  -  linear,  C  -  circular) 

1 

!  Mode  of  calculating  the  radiation  pattern 

Figure  6.22  -  Input  file  cimpa.in 


First  the  CyMPA  routine  was  run  to  calculate  the  current  distribution  at  each  patch  and  the  field 
distribution  in  each  aperture  in  the  array,  as  well  as  mutual  coupling  and  input  impedance,  in  the 
frequency  band  of  interest  4.5-6.5  GHz.  In  the  output  file  “Sij.daf'  the  S-parameters  are  written 
as  a  function  of  frequency: 


#  GHz  S  DB  R  50 
4.5000  -2.85  -144.19 

-24.33  -6.86 

-21.89  -63.58 

-32.05  -123.28 
4.6000  -3.99  -159.57 

-23.29  -31.30 

-22.63  -76.29 

-30.13  -147.50 
-5.62  -176.97 
-22.65  -57.39 

-23.29  -84.31 

-28.48  -172.52 


-24. 

.33 

-6. 

.86 

-21. 

.89 

-2. 

.85 

-144. 

.19 

-32. 

.05 

-32. 

.05 

-123. 

.28 

-2. 

.85 

-21. 

.89 

-63. 

.58 

-24. 

.33 

-23. 

.29 

-31. 

.30 

-22. 

.63 

-3. 

.99 

-159. 

.57 

-30. 

.13 

-30. 

.13 

-147, 

.50 

-3. 

.99 

-22. 

.63 

-76. 

.29 

-23. 

.29 

-22, 

.65 

-57. 

.39 

-23. 

.29 

-5. 

.62 

-176. 

.97 

-28. 

.48 

-28. 

.48 

-172. 

.52 

-5, 

.62 

-23. 

.29 

-84. 

.31 

-22. 

.65 

-63, 

.58 

-32 

.05 

-123. 

.28 

-123, 

.28 

-21, 

.89 

-63. 

.58 

-144, 

.19 

-24, 

.33 

-6. 

.86 

-6, 

.86 

-2 

.85 

-144. 

.19 

-76. 

.29 

-30, 

.13 

-147, 

.50 

-147. 

.50 

-22 

.63 

-76. 

.29 

-159. 

.57 

-23, 

.29 

-31. 

.30 

-31. 

.30 

-3, 

.99 

-159. 

.57 

-84. 

.31 

-28. 

.48 

-172. 

.52 

-172. 

.52 

-23. 

.29 

-84. 

.31 

-176. 

.97 

-22 

.65 

-57, 

.39 

-57. 

.39 

-5. 

.62 

-176. 

,97 

4.7000 
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The  comparison  of  calculated  and  measured  Sn  and  the  S21  parameters  is  shown  in  Figs.  6.23 
and  6.24.  Notice  that  the  peaks  and  drops  in  measured  Su  and  S21  curves  for  frequencies  larger 
than  6  GHz  are  due  to  the  resonances  in  the  feeding  structure,  i.e.  they  are  due  to  the  finite  length 
of  the  antenna  in  the  axial  direction. 


Figure  6.23.  Calculated  and  measured  magnitude  ofSu. 


Figure  6.24.  Calculated  and  measured  magnitude  0/821  for  patches  in  azimuthal  direction. 
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The  CyMPApat  routine  was  run  to  calculate  the  radiation  pattern.  Since  we  wanted  radiation 
pattern  (with  mutual  coupling  taken  into  account)  at  5.2  GHz,  we  changed  the  first  line  of  the 
input  file  “cimpa.in”  to  run  this  case: 

5.2  6.5  11  R !  Frequency  (GHz) 

The  comparison  between  calculated  and  measured  radiation  pattern  is  given  in  Fig.  6.25. 


azimuthal  angle  [deg] 


Figure  6.25.  Radiation  pattern  of  2x2  array  of  aperture-coupled  patches  in  azimuthal  plane 
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7  SMiPAWin  and  CympaWin 
USER’S  MANUAL 


graphical  user  interface  for 

SMiPA 

SMiPApat 


CyMPA 

CyMPApat 
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SMiPAWin  and  CympaWin 
USER’S  MANUAL 


graphical  user  interface  for 

SMiPA 

SMiPApat 

CyMPA 

CyMPApat 


Introduction 

The  SMiPA,  CyMPA,  SMiPApat  and  CyMPApat  programs  have  a  user  interface  through  text 
input  and  output  files,  for  both  input  and  output.  For  the  convenience  of  their  users  working  on 
PCs  using  Microsoft  Windows®  the  SMiPAWin  and  CympaWin  software  packages  are 
developed.  SMiPAWin  package  includes  the  SMiPA,  SMiPApat  and  SMiPAWin  -  a  Windows 
based  graphical  user  interface.  Simmilarly,  CympaWin  software  package  includes  CyMPA, 
CyMPApat  and  CympaWin  -  Windows  based  graphical  user  interface.  SMiPAWin  and 
CympaWin  run  on  IBM  compatible  PC  with  at  least  386  processor,  16  Mb  RAM  and  Windows 
95/98  or  later  operation  system  installed.  SMiPAWin  and  CympaWin  are  very  simmilar  fi'om  the 
users  point  of  view,  enabling  identical  features  for  spherical  (SMiPAWin)  and  cylindrical 
(CympaWin)  curved  microstrip  structures.  Therefore,  features  of  these  two  software  packages 
are  presented  in  parallel,  mentioning  both  of  them  at  each  described  feature.  For  the  convenience 
the  SMiPAWin  is  and  the  corresponding  CympaWin  equivalent  is  mentioned  in  brackets. 

SMiPAWin  (CympaWin)  enables  the  following: 

-  easy,  spreadsheet-like  creating  and  editing  of  the  input  file  smipa .  in  (cympa .  in  ) 

-  launching  of  the  simmulation  directly:  SMiPAWin  (CympaWin)  executes  SMiPA  (CyMPA) 
and  SmiPApat  (CyMPApat)  from  MS  -  DOS  prompt  automatically 

-  display  of  calculated  S-parameters,  input  impedance  or  radiation  pattern  using  data  from 
output  files 
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These  features  enable  the  user  to  create  input  files  for  SMiPA  (CyMPA)  and  SmiPApat 
(CyMPApat)  with  a  click  of  a  mouse,  without  need  for  careful  text  editing  within  smipa .  in 
(cympa .  in)  file. 


The  Installation  of  the  Software 

Execute  setup .  exe  from  SMiPA  (CyMPA)  installation  disk  or  CD  ROM.  During  instalation 
you  can  decide  in  which  folder  you  want  to  have  your  SMiPAWin  (CympaWin),  SMiPA 
(CyMPA)  and  SmiPApat  (CyMPApat)  installed  (e.g.  in  folder  c:\Program  Files\Smipa 
or  c :  \Program  Files\Cympa). 

Now  you  can  start  the  SMiPAWin  (CympaWin),  entering  your  first  antenna  configuration  and 
running  calculations. 


Using  SMiPAWin  (CympaWin) 

The  flow  chart  for  SMiPAWin  and  CympaWin  shows  interactions  between  all  identities 
involved;  user,  SMiPAWin  (CympaWin),  SMiPA  (CyMPA)  and  SmiPApat  (CyMPApat) 
programs.  SMiPAWin  (CympaWin)  enables  saving  the  configuration  in  arbitrary  folder  we  call 
design  folder,  under  the  project  name  (here  pro  ject_n£une  will  be  used)  selected  through  the 
standard  dialog  box.  All  output  files  will  be  named  project_name  plus  corresponding 
extension.  The  configuration  file  extension  is  'sin'  ('cin'  for  CympaWin),  so  complete  name  of 
the  configuration  file  is  pro  ject_naine .  sin  (pro  ject_name .  cin). 

Here  are  listed  some  hints  that  give  more  in  depth  information  about  using  and  operation  of 
SMiPAWin  (CympaWin): 

•  The  SMiPA  (CyMPA)  and  SMiPApat  (CyMPApat)  use  input  file  smipa .  in  (cympa .  in) 
where  all  configuration  information  is  stored,  smipa .  in  (cympa .  in)  file  is  created  by 
SMiPAWin  (CympaWin)  as  a  copy  of  the  pro  ject_name .  sin  (pro  ject_name .  cin) 
file.  After  calculation,  output  files  with  calculation  results  calculated  by  SMiPA/SMiPApat 
(CyMPA/CyMPApat)  occur  in  design  folder. 

•  Rigorous  calculation  is  performed  by  execution  of  SMiPA  (CyMPA).  Thus,  input  impedance 
and  S  parameters  are  obtained. 

•  SMiPApat  (CyMPApat)  can  calculate  radiation  pattern  using  rigorous  models  {Gain  pattern, 
Pattern  normalized  to  reference  direction,  Element  pattern)  or  simple  models  {Simple  model 
-  uniform  excitation.  Simple  model  -  file  for  patch  excitation,  Simple  model  -  maximum 
predefined phases  corrected).  CAUTION:  For  rigorous  models  SMiPA  (CyMPA)  must  be 
run  first! 

•  It  is  not  possible  to  run  calculations  if  the  configuration  was  changed  without  saving  it  first. 
The  application  title  bar  indicates  the  file  name,  and  '*'  is  appended  if  file  was  changed  but 
not  saved. 
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•  It  is  possible  to  view  the  previous  results  without  running  the  calculations  again.  When  the 
configuration  name  is  selected,  SMiPAWin  (CympaWin)  will  enable  View  submenu  items  if 
previous  output  results  of  the  configuration  calculation  were  found  in  the  selected  design 
folder. 

•  If  new  results  are  needed,  configuration  should  be  edited,  saved,  and  SMiPA  (CyMPA) 
and/or  SMiPApat  (CyMPApat)  programs  should  be  executed  according  to  the  SMiPAWin 
(CympaWin)  flow  chart.  Then,  new  results  can  be  viewed  under  the  View  menu.  Graphs  of 
calculation  results  will  be  shown  after  calculation  execution  when  selected  in  View  menu. 


Input  and  output  files 

SMiPA  (CyMPA)  and  SMiPApat  (CyMPApat)  both  use  fixed  file  names  and  will  be  saved  after 
calculation  into  the  folder  we  call  design  folder  as  follows: 

temp2 .  t  output  of  SMiPA  (CyMPA)  containing  the  patch  mutual  coupling  and  input 
impedance  information,  SMiPAWin  (CympaWin)  copies  it  automatically  as 
pro ject_name . si j 

pat-az .  dat  output  of  SMiPApat  (CyMPApat)  containing  field  values  in  the  azimuthal  plane, 
SMiPAWin  (CympaWin)  copies  it  automatically  as  pro  ject_naine  .paz 

pat-el .  dat  output  of  SMiPApat  (CyMPApat)  containing  field  values  in  the  elevation  plane, 
SMiPAWin  (CympaWin)  copies  it  automatically  as  pro  ject_name  .pel 

These  three  files  contain  sufficient  information  about  the  simulation  run. 

Also  the  various  temporaiy  files  will  be  generated  during  the  program  run  within  the  design 
folder.  Some  of  them  are  used  for  status  information  and  synchronization  between  the  Windows® 
and  DOS  application  and  are  deleted  before  or  after  the  program  execution. 


The  user  interface  properties 

•  SMiPAWin  (CympaWin)  is  designed  as  the  standard  Windows  SDI  (Single  Document 
Interface).  It  maintains  the  look  and  feel  of  standard  Windows  applications  to  minimize  the 
learning  curve. 

•  In  order  to  minimize  typing  errors  SMiPAWin  (CympaWin)  forces  user  to  enter  correct  type 
of  parameters,  and  forbids  saving  configuration  with  empty  values  in  forms. 

•  The  SMiPAWin  (CympaWin)  is  menu  oriented  with  submenu  items  logically  grouped,  and 
disabled  ('grayed')  if  operation  is  not  possible  within  the  current  context.  Standard  dialog 
boxes  are  used  for  the  file  operations. 

•  The  Toolbars  exist  which  duplicate  some  more  frequent  operations. 
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generated  in  design  folder: 
pro j  ect_f ile . out 
pro ject^file.  pirn 
project_file. si j 

content  displayed: 
pro ject_f ile . out 

Port  impedance  (.niml  and  S- 
parameter  (.sii)  graphs  enabled! 

jiBle''Ydiryiew' 

rola^lol  #1 1  l(i 


generated  in  design  folder: 
pro ject_f ile . out 
project_file.paz 
pro  j  ec t_f i le . pel 

content  displayed: 
smipapat .out 
(cympapat . out) 

Azimuth  (.paz)  and  elevation  (.pel  > 
radiation  pattern  graphs  enabled! 


SMiPAWin 
I  (CympaWin) 


execution  of 
SMiPA 
CyMPA 


rigorous 

pattern 


execution  of 
SMiPApat 
(CyMPApa) 


approximate  pate(n 

calculation  | 

(Simple  models  only:  | 

Edit  Configuration  \ 


i  jile  "Edit  yiew  flun  "Help 

[  DtelHl  «l  !  di'l 


SMiPAWin  (CympaWin)  Flow  Chart:  interactions  between  entities  involved 
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Some  SMiPAWin  (CympaWin)  windows 

In  this  section  some  typical  windows  that  occur  during  the  use  of  SMiPAWin  (CympWin)  are 

shown  and  explained.  Also,  'Run'  and  'View'  menues  are  described. 

Window  1 

File  New  or  File  ->  Open 

•  This  should  be  the  first  user  action,  since  most  of  the  other  operations  are  not  possible  until 
the  configuration  name  is  known. 

•  Window  1  shows,  that  the  configuration  file  'Testcin'  is  already  loaded  and  some  of  the 
toolbars  are  enabled. 

Window  2 

Configuration  creating,  editing,  and  viewing 

•  Realized  with  the  Page  Control  Dialog  Box  (Style  Sheets)  as  shown  in  Window  2 

•  Configuration  parameters  are  grouped  logically  and  placed  at  the  same  card 

Window  3 

Error  message  box 

•  Any  parameter  not  entered  will  cause  the  error  message  box  to  pop  up  when  user  tries  to 
leave  Configuration  dialog  box  using  ‘Apply’  button 

•  It  is  possible  to  leave  the  control  without  filling  it  completly,  by  using  ‘Cancel’  button 

Window  4 

•  Confirmation  will  be  requested  if  user  tries  to  leave  the  program  without  saving 
configuration  which  was  edited/changed. 

•  The  application  does  not  ask  for  confirmation  when  ‘Cancel’  button  is  used,  since  it  is 
always  possible  to  go  back  to  editing. 

Window  5 

•  Tables  resize  automatically  when  parameters  that  determine  the  number  of  rows  change. 

•  Scroll  bars  are  shown  when  table  does  not  fit  the  page  control  size. 

Run  Menue 

•  Run  menu  contains  two  submenu  items:  Execute  SMiPA  {Execute 
Cympa)  and  Execute  SMiPApat  {Execute  CympaPat). 

•  SMiPA  (Cympa)  and  SMiPApat  (CympaPat)  applications  are 
executed  as  minimized  DOS  boxes. 

•  User  can  maximize  them  and  watch  the  program  output  as  within  the  command  line 
environment. 

View  Menue 

•  Configuration  submenue  is  enabled  first  when  configuration  is 
created  by  File  ->  New  or  reload  using  File  ->  Open. 

•  S-parameters  graph  and  Port  impedances  graph  are  enabled  after 
execution  of  SMiPA  (CyMPA) 

•  Radiation  pattern  greaph  is  enabled  after  execution  od 
SMiPApat  (CyMPApat) 
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Troubleshooting 

If  SMiPA  (CyMPA)  or  SMiPApat  (CyMPApat)  would  not  start,  try  running  them  directly  from 
the  command  line  within  the  folder  where  they  are  installed  (e.g.  in  folder  c:\Program 
Files\Smipa  or  c :  \ Program  Files\Cympa).  This  folder  must  contain  the  smipa .  in 
(cympa .  in)  file  that  was  generated  by  SMiPAWin  (CympaWin).  Therefore  you  must  copy  it 
from  design  folder  before  running  desired  program.  The  output  files  will  in  this  case  appear  in 
the  same  folder  and  will  be  named  corresponding  to  names  in  section  ‘input  and  output  files’. 
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Conclusions 


The  outcome  of  the  proposed  two-phase,  24-month  effort  was  developing  a  model  for  analyzing 
single-curved  and  double-curved  conformal  microstrip  antennas.  Furthermore,  it  was  planned  to 
develop  two  computer  programs:  (a)  for  analyzing  patch  arrays  printed  on  spherical  structures, 
and  (b)  for  analyzing  cylindrical  aperture-coupled  patch  arrays. 

We  have  developed  the  program  SMiPA  for  analyzing  microstrip  patch  arrays  printed  on 
spherical  structures.  It  is  assumed  that  the  patches  have  either  rectangular  or  circular  shape. 
Three  types  of  feeding  structures  are  considered:  coaxial  and  microstrip  transmission  lines  and 
aperture  coupling.  The  program  solves  the  electric  field  integral  equation,  and  the  moment 
method  is  used  for  solving  the  integral  equation.  The  analysis  of  patch  antennas  printed  on  (or 
embedded  in)  multilayer  spherical  structures  is  obtained  by  implementing  a  routine  for 
calculating  Green's  functions  of  multilayer  spherical  structures.  The  numerical  evaluation  of 
matrix  elements  used  in  moment  method  is  made  with  special  care.  In  particular,  we  defined  the 
normalized  Legendre  polynomials  and  vector-Legendre  transformation  by  which  evaluation  of 
ratios  of  very  large  numbers  was  avoided,  and  numerical  treatment  of  spherical  Bessel  and 
Hankel  functions  was  carefully  made. 

The  program  SMiPA  calculates  the  following  parameters  of  the  cylindrical  patch  arrays: 

•  input  impedance  at  each  input  port  in  the  array 

•  mutual  coupling  between  each  two  patches  in  the  array 

•  radiation  pattern  of  the  array  when  all  mutual  couplings  are  taken  into  account 

•  radiation  pattern  of  the  array  without  taking  mutual  coupling  into  account  (fast 
calculations  of  the  radiation  pattern  are  needed  for  making  first  design  of  the  array). 

Furthermore,  we  have  developed  the  upgrade  of  the  program  CyMPA  that  calculates  aperture- 
coupled  patch  antennas  on  cylindrical  structures.  Program  CyMPA  is  the  outcome  of  the  project 
F61 775-99- WE040,  “CyMPA  -  Program  for  Analyzing  Microstrip  Patch  Arrays  on  Circular- 
Cylindrical  Structures,”  sponsored  by  BOARD.  Old  version  of  the  program  enables  the  analysis 
of  cylindrical  patch  arrays  fed  by  a  coaxial  or  microstrip  transmission  line.  With  the  upgraded 
program  it  is  possible  to  fully  analyze  cylindrical  patch  arrays  with  aperture  feeding,  which  is 
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important  since  aperture  type  of  feeding  enables  one  to  separate  radiating  and  feeding  part  of  the 
antenna,  thus  enabling  independent  design  both  of  antenna  parts.  The  upgraded  program  fully 
analyzes  the  array,  i.e.  radiation  pattern,  input  impedance  at  each  input  port,  and  mutual  coupling 
between  each  two  elements  is  calculated. 

The  programs  are  tested  by  comparing  the  calculated  results  with  measurements  and  with 
calculated  results  obtained  by  other  methods.  We  have  additionally  developed  an  experimental 
model  for  testing  the  program  for  analyzing  spherical  arrays  since  there  are  only  a  few 
experimental  results  in  literature  in  connection  with  spherical  patch  antennas.  In  all  tested  cases 
there  is  a  good  agreement  between  calculated  results  and  test  data. 
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